Improving the sensitivity to gravitational-wave sources by modifying the input-output 

optics of advanced interferometers 
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To improve the sensitivity of laser-interferometer gravitational-wave (GW) detectors, experimen- 
tal techniques of generating squeezed vacuum in the GW frequency band are being developed. The 
squeezed vacuum generated from non-linear optics have constant squeeze angle and squeeze factor, 
while optimal use of squeezing usually requires frequency dependent (FD) squeeze angle and/or ho- 
modyne detection phase. This frequency dependence can be realized by filtering the input squeezed 
vacuum or the output light through detuned Fabry-Perot cavities. In this paper, we study FD input- 
output schemes for signal-recycling interferometers, the baseline design of Advanced LIGO and the 
currently operational configuration of GEO 600. Complementary to a recent proposal by Harms 
et al. to use FD input squeezing and ordinary homodyne detection, we explore a scheme which 
uses ordinary squeezed vacuum, but FD readout. Both schemes, which are sub-optimal among all 
possible input-output schemes, provide a global noise suppression by the power squeeze factor. At 
high frequencies, the two schemes are equivalent, while at low frequencies the scheme studied in 
this paper gives better performance than the Harms et al. scheme, and is nearly fully optimal. We 
then study the sensitivity improvement achievable by these schemes in Advanced LIGO era (with 
30-m filter cavities and current estimates of filter- mirror losses and thermal noise), for neutron star 
binary inspirals, for low-mass X-ray binaries and known radio pulsars. Optical losses are shown 
to be a major obstacle for the actual implementation of these techniques in Advanced LIGO. On 
time scales of third-generation interferometers, like EURO/LIGO-III (~ 2012), with kilometer-scale 
filter cavities and/or mirrors with lower losses, a signal-recycling interferometer with the FD readout 
scheme explored in this paper can have performances comparable to existing proposals. 

PACS numbers: 04.80.Nn, 03.65.ta, 42.50.Dv, 95.55.Ym 



I. INTRODUCTION 

The first generation of kilometer-scale, ground-based 
laser-interferometer gravitational-wave (GW) detectors 
(interferometers for short), located in the United States 
(LIGO 0,0), Europe (VIRGO and GEO 600 00) 
and Japan (TAMA 300 0), have begun their search for 
gravitational radiation and have yielded first scientific 
results 0, H> The development of interferometers 
of the second generation, such as Advanced LIGO (to 
be operative around 2008 [10]L and future generations 
(such as EURO and LIGO-III), is underway. In this pa- 
per we explore the possibility of improving the sensitiv- 
ity of signal-recycling (SR) interferometers 0, 0] , the 
baseline design of Advanced LIGO ^} an d the current 
optical configuration of GEO 600 , when squeezed vac- 
uum is injected into the antisymmetric port (the "input 
port", as we shall refer to it in this paper 1 ). 

In the early 1980s, building on works of Caves [l3[ . 
Unruh Q proposed the first design of a squeezed- 
input interferometer, which can beat the free-mass Stan- 
dard Quantum Limit (SQL) 0|. Other theoretical 
studies of input squeezing followed 16] . If generated 



from non-linear optics, squeezed vacuum will have fre- 
quency independent squeeze angle and squeeze factor in 
the GW frequency band 0]. The above theoretical 
works, as well as past experiments employing squeezed 
vacuum to enhance interferometer performances |18| . 
all assume frequency independent squeezing. In the 
1990s, Vyatchanin, Matsko and Zubova 0] realized that 
the sensitivity of GW interferometers can also be im- 
proved, beating the SQL, by measuring an optimal out- 
put quadrature, which is usually frequency dependent. 
Later, Kimble, Levin, Matsko, Thorne and Vyatchanin 
(KLMTV) |2Cj made a comprehensive, unified theoret- 
ical study of improving the sensitivity of conventional 
interferometers 2 by injecting squeezed vacuum into the 
input port and/or performing frequency dependent (FD) 
homodyne detection at the output port. They showed 
that, for conventional interferometers, in order to obtain 
a noise suppression proportional to the power squeeze fac- 
tor at all frequencies (optimal input squeezed vacuum), 
either the squeezed vacuum must have a FD squeeze an- 
gle {squeezed-input interferometer), or FD homodyne de- 
tection has to be applied at the output port (squeezed- 
variational interferometer). [Of course, combinations of 



1 This is the same port from which the GW signal light exits, but 
here the squeezed vacuum propagates into the interferometer, 
instead of coming out of it. 



2 By conventional interferometer we mean a Michelson interfer- 
ometer without SR cavity, or with SR cavity on resonance or 
antiresonance with the laser frequency. 
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those optical configurations can also be used, but they 
will be experimentally more challenging.] KLMTV pro- 
posed a practical way of implementing FD homodyne 
detection, as well as converting squeezed vacuum with 
constant squeeze angle into squeezed vacuum with FD 
squeeze angle, by filtering the output light or input 
squeezed vacuum through two detuned (with respect to 
the laser carrier frequency) FP cavities (KLMTV filters). 
KLMTV constructed the explicit filter parameters that 
provide the desired frequency dependence for squeezed- 
input and squeezed-variational interferometers, showing 
that the latter provides a better ideal performance than 
the former, but is more susceptible to optical losses. 
Purdue and Chen (PC) studied the KLMTV filters fur- 
ther, and worked out the most general FD squeeze angle 
and homodyne phase that a sequence of filters can pro- 
vide [m. 

Experimental programs on generating squeezed vac- 
uum in the GW frequency band and injecting it into GW 
interferometers have already started in several groups, 
for example, at the Australian National University [22J, 
at the Massachusetts Institute of Technology, USA [23J , 
and at the Albert Einstein Institut in Hannover, Ger- 
many [2^|. Their goal in the next several years is to 
inject squeezed vacuum with ~10 dB squeeze factor (as a 
net result after optical losses) into an interferometer. It 
is very likely that the squeezed vacuum they obtain has 
a constant squeeze angle in the GW frequency band. 

This paper contains three relatively independent parts. 
In the first part, we generalize the work by KLMTV 
on FD input-output optics to SR interferometers. Re- 
cently, Harms et al. [25j applied the KLMTV squeezed- 
input scheme, combining FD input squeezed vacuum 
with ordinary homodyne detection to SR interferometers, 
achieving a global noise suppression equal to the power 
squeeze factor. Harms et al. also showed that their FD 
squeezed-input scheme is only sub-optimal; the fully op- 
timal scheme, however, requires complicated frequency 
dependence in both input squeeze angle and homodyne 
phase, and cannot be achieved by KLMTV filters. Com- 
plementary to Harms et al.'s scheme, we explore here 
the scheme which combines ordinary input squeezed vac- 
uum with FD homodyne detection (henceforth, the BC 
scheme) . This scheme, which can be thought of as a gen- 
eralization of the KLMTV squeezed-variational scheme, 
can also provide a global noise suppression by the power 
squeeze factor. In addition, at high frequencies (above 
~ 200Hz for typical Advanced LIGO configurations), it 
is equivalent to the Harms et al. scheme, while at low 
frequencies (below ~ 200 Hz) it is to a very good ap- 
proximation fully optimal, and thus provides a better 
sensitivity than the Harms et al. scheme. 

In the second part of this paper we apply these FD 
input-output schemes to Advanced LIGO (2008), assum- 
ing that the generation and injection of squeezed vac- 
uum might have already (or partially) become available 
at that time scale. The major obstacle in using FD input- 
output techniques in the facilities of Advanced LIGO is 



the constraint that the filter cavities cannot be longer 
than ~30 meters — the shorter the filter cavities, the 
larger the optical losses. In our analyses we assume that 
filter losses dominate over internal interferometer losses, 
and comment only briefly on the effects of the latter. To 
quantify the improvement in sensitivity due to FD tech- 
niques, we consider three classes of astrophysical sources: 
neutron-star (NS) binary inspirals, low-mass X-ray bi- 
naries (LMXBs) and known radio pulsars. In addition 
to the ideal quantum noise and filter optical losses, we 
also take into account current estimates of thermal and 
seismic noises. [We note that GW interferometers can 
already take advantage of input squeezed vacuum even 
without using FD techniques, if the interferometer pa- 
rameters are carefully optimized. For example, an in- 
teresting optical configuration without FD input-output 
optics has been explored by Corbitt and Mavalvala |23| . 
providing good sensitivities at high frequencies.] 

In the third part of the paper, we apply our FD read- 
out scheme to third-generation interferometers, such as 
EURO/LIGO-III, which are scheduled to be operative 
around 2012. We assume that on this time scale, due to 
the implementation of cryogenic techniques and the use 
of kilometer- length KLMTV filters, thermal noise will 
be negligible and loss effects will be rather low. Third- 
generation interferometers will have to beat the SQL sig- 
nificantly. We compare the performance of SR interfer- 
ometers with our FD readout scheme with those of other 
existing SQL-beating proposals, such as the KLMTV 
squeezed-variational interferometer |2(j and the speed- 
meter interferometers [2(| [23). We also investigate the 
accuracy of short-arm and short-filter approximations 
used in describing GW interferometers and KLMTV fil- 
ters. More dramatic ideas to circumvent the SQL in 
GW interferometers exist, for example, the intra-cavity 
schemes of Braginsky, Gorodetsky and Khalili |28| . and 
the feed-back control scheme of Courty, Heidmann and 
Pinard [2^. Since thorough analyses of these schemes 
tuned to GW interferometers has not been available yet, 
in this paper we do not compare the performances. 

Readers with particular interests in the astrophysical 
consequences of using input squeezed vacuum and FD 
schemes could go directly to the second part of the pa- 
per [Sec. 0, in which an in-depth understanding of the 
optics is not required. The paper is organized as fol- 
lows. In Sec. |Hj we write the input-output relation of 
a non-squeezed SR interferometer in terms of the intrin- 
sic FD rotation angle and ponderomotive squeeze factor. 
In Sec. IIHI we review the KLMTV filters, including the 
effects of optical losses. In Sec. IIVI we study FD input- 
output schemes for SR interferometers. More specifically, 
in Sec. IIV Al we write the general input-output relation 
of SR interferometers with FD input-output optics; in 
Sec. IIVBI we study all sub-optimal schemes that allow 
global noise suppression, proposing the BC scheme; in 
Sec. ESI we study the regime with low ponderomotive 
squeezing (high frequency band of Advanced LIGO), and 
show the equivalence between the BC and the Harms 
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Quantity 


Symbol and Value 


Laser frequency 


u>o = 1.8 x lO^sec" 1 


GW sideband frequency 


Q, 


Arm-cavity length 


L (4 km for LIGO facilities) 


Mirror mass 


m (40 kg for Advanced LIGO) 


Input test-mass (ITM) power transmissivity (LIGO only) 


T 


Arm-cavity circulating power 


I c (840 kW for Advanced LIGO) 


Light power at the beamsplitter 


Io 


SR optical resonant (sideband) frequency 


— A 


SR bandwidth 


c 


Homodyne detection phase 


c 


Input squeeze factor 


r 


Input squeeze angle 


a 



TABLE I: Parameters of the SR interferometer and input squeezed vacuum. 



et al. schemes; in Sec. II V Dl we study the fully optimal 
scheme, showing that at low frequencies the BC scheme 
is a good approximation to it. In Sec. we investi- 
gate the improvement in sensitivity to GWs from various 
astrophysical sources. In Sec. IVII we compare the BC 
scheme with other proposals for third-generation inter- 
ferometers, and study the effect of filter lengths in FD 
readout schemes. Finally, Sec. IVIll summarizes our con- 
clusions. 



II. QUADRATURE ROTATION AND 
PONDEROMOTIVE SQUEEZING IN SIGNAL 
RECYCLED INTERFEROMETERS 



A summary of the various parameters of SR interfer- 
ometers, such as Advanced LIGO and GEO 600, is given 
in Table [H The input-output relation for the quadra- 
ture fields in signal recycled interferometers reads [see, 
e.g., Eq. (24) in Ref. |30|. with superscript (1) and tilde 



dropped] 



1 

M 



Cn C12 

C21 C22 



a 1 

«2 



D 2 



hsQL 



where we define 



M= [A 2 - (ft + ief] n 2 -\t c , 



and 



C11 = C22 
C12 
C21 



n 2 (fl 2 -\ 2 +e 2 ) + \i c 
-2e\fl 2 , 
2e\fl 2 ~ 2ei c , 



(i) 



(2) 



(3) 
(4) 
(5) 



(G) 



The parameters A and e are related to the real and imag- 
inary parts of the free 3 optical resonant frequency Wf rcc 
of the SR interferometer by [3(| : 



u4 r B ee = ujq - A - ie . 



where loq is the laser frequency. The parameter t c is 
defined by 



8uj I c 
mL c 

(2tt x 100 Hz) 3 



840 kW 



40 kg 



4 km 
L 



,(8) 



where c is the speed of light, m is the mirror mass, L is 
the arm length, and I c is the circulating optical power in 
the arm cavity, which in turn depends on the power at 
beamsplitter, /qi by 4 



Ic = y Jo • 



(9) 



The quantity l c must be on the order of ^q W if we want 
the opto-mechanical coupling to modify the detuned (as- 
suming A ~ ficw) interferometer's sensitivity in the GW 
frequency band [see Eqs. ©~©]- In addition, we have 
denoted by 



^SQL 



87; 



mn 2 L 2 



(10) 



the free-mass SQL for the gravitational strain /i(O). 5 

It is important to note that the input-output relation, 
as given by Eqs. 0~©j has been obtained at the lead- 

1 /3 

ing order in QL/c (as well as in eL/c, XL/c and l c L/c). 



3 Here "free" means that the mirrors are all fixed at their equilib- 
rium positions. 

4 This only applies to LIGO; for GEO 600 we have I c = I /2. 

5 Note that the definitions of t c and /isql, written in terms of m, 



L and I c , 
GEO 600. 



differs by numerical factors in Advanced LIGO and 
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This approximation is called the "short-arm" approxi- 
mation, since it assumes that the arm length be much 
smaller than the gravitational wavelength. The short- 
arm approximation simplifies dramatically the form of 
the input-output relation, as well as the design of opti- 
mal KLMTV filters (as we shall see later in this paper). 

In the following sections we first write the input-output 
bi.2-o-i.2 relation [the first term inside the parenthesis on 
the RHS of Eq. QJ] in terms of an intrinsic squeeze factor 
q and an intrinsic rotation angle p. We then study how 
the output quadratures depend on the signal [the second 
term inside the parenthesis on the RHS of Eq. QJ] , ob- 
taining the quadrature Cmax at which the signal strength 
is maximal. Finally, we give the noise spectrum, and 
express it in terms of q, ip and Cmax- 



9 must be the rotation angle of M; otherwise, in order 
that R(0)S(g, p) = R(9')S(q, p'), one must impose that 



(17) 



where m and n are integers, and m — n must be an even 
number. The uniqueness of the decomposition i|12|) as- 
sures that the quantities 9, q and <p have unambiguous 
physical meaning. 

By comparing Eq. (|12[1 with Eqs. 10 -0, we can ex- 
press the angles 9, <p and the factor q in terms of SR 
parameters. Using the identity (JTJJ, we get the squeeze 
factor q: 



cosh2g = 1 



M 



(18) 



A. Rotation of noise quadratures and 
ponderomotive squeezing 

As is evident from Eqs. m detuned SR inter- 

ferometers (with A 7^ 0), the input-output relation is fre- 
quency dependent. For high-power interferometers like 
Advanced LIGO, the matrix dj contains both a FD ro- 
tation and a FD (ponderomotive) squeezing. Let us work 
these quantities out explicitly. The quantum transfer ma- 
trix (with an overall phase factor removed) 



i, j" = 1,2 (detA% = l) 



(11) 



is a matrix with real elements and with determinant equal 
to one |3lj| . So, it can always be written in the form 



M = R(9)S(p,q). 



(12) 



as a product of a rotation operator R(#) and a squeezing 
operator S(ip, q), defined by 



R(0) = 



cos 9 



- sin 9 
cos 6 



and 



s(«) 



R(p)S(q)K(-p), 



q>0. 



(13) 

(14) 
(15) 



Here 9 is the rotation angle, q is the squeeze factor, and 
ip is the squeeze angle. These quantities can in general 
be frequency dependent. It can be easily shown that the 
decomposition l|12|) is also unique, unless M. is a pure 
rotation (with q = 0, 9 being the rotation angle, and ip 
being arbitrary). From Eqs. (|12fl - l|14[) . we have 



Ti(MM^) = 2cosh2g, 



(16) 



which determines q uniquely once M is given (recall that 
q > 0). If q is zero, S(y>, q) reduces to the identity matrix 
regardless of the value of <p [see Eqs. I|14fl and l|15f) ]. and 



Since for SR interferometers, C\\ — C22, we must impose 
cos(2<,9 + 9) = 0, or 6 = tt/2 — 2ip. This casts M.^ into 
the following form: 

M 

= R(n/2-<p)S(q)R(-<p) 

I sin 2p cosh q — sinhg — cos2</? coshg \ 
I sinhg — cos 2(^9 coshg sin 2p cosh q J 

(19) 



Thus we have 



tan 2p = — 



Mn + M22 



sinhq 



M12 + M21 
l c X + Q 2 (e 2 - A 2 + ft 2 ) 

e(2A^ 2 - l c ) 
M21 - M12 _ eic_ 
2 ~ \M\ ' 



(20) 
(21) 



Note that Eq. (2T) agrees with Eq. I|18|) . 

In absence of SR mirror, or when the SR cavity is either 
resonant or antiresonant with the carrier, we have A = 0, 
and the above equations reduce to the known expressions 
for conventional interferometers 201 



sinh q 



K 
~2 



with 



K = 



tan 2<p 



2eir 



2 

"jC' 



n 2 (n 2 + e 2 ) 



(22) 



(23) 



K, being the coupling constant defined by KLMTV in 
Eq. (18). The smhq in Eqs. J2U and (|2"2|) is proportional 
to l c , which is in turn proportional to the circulating 
power I c and inversely proportional to the mirror mass 
m. This means the squeezing arises from the well-known 
ponderomotive effect [32j. 

In Fig. ^ we pl°t 2 sinh q (left panel) and 2p> (right 
panel) as functions of frequency, for two typical SR con- 
figurations and a conventional-interferometer configura- 
tion. [Note that 2 sinh q ~ 2q as q <C 1 and 2 sinh q ~ e q 
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FIG. 1: We plot 2 sinh g (left panel) and 2ip (right panel) for two typical SR optical configurations, the broadband with 
A = 2tt x 234.07 Hz, e = 2tt x 70.36 Hz, (dashed line) and narrowband with A = 2tt x 601.43 Hz, e = 2tt x 25 Hz, (dot-dashed 
line). We also plot the curves for the conventional interferometer with A = 0, e = 2-7T x 93.75 Hz (continuous line). In all cases 
we fix I c = 840 kW and m = 40 kg. 



as g > 1.] As we can see from the plots, both 2 sinh g 
and 2ip are frequency dependent. 

Let us focus on the detuned configurations (dash and 
dash-dot curves in Fig.^l. The squeeze factor decreases 
at high frequencies. This can be easily understood from 
Eq. J2U, where we see that sinhg (hence q) decreases 
when f2 increases, because M, Eq. (J2J, is a polynomial 
in fi, so it grows indefinitely as f2 tends to infinity. [The 
factor r2 2 in Eq. © can be traced back to the response 
of a free mirror to an external force, which decreases as 
l/f2 2 ; while the factor [A 2 — (fl + ie) 2 ] increases at high 
frequencies because the storage time (1/e) of the inter- 
ferometer becomes much longer than the GW period.] 
Using Eq. © and the fact that \M\ > 3(M), we obtain 

\M\>2efl 3 . (24) 

Combining Eq. with Eqs. lj2T)l and (JSJ, we have 

sinh q = - — % < — ^7 



2 x 10" 



2tt x 300 Hz 

n 



40kg\ /4km 



m J \ L 



840 kW, 
(25) 



cal spring" resonance 6 (left peaks), and mildly near the 
optical resonance (right peaks). Those resonant features 
in q are caused by local minima of M around the two 
resonant frequencies; the optical resonance provides less 
squeezing since squeezing is already suppressed at such 
high frequencies [see Eq. (|2"5)) ]. The squeeze factor tends 
to a nonzero constant for Q much lower than the resonant 
frequencies. By taking the limit of Eq. (|21|l when Q — > 
we obtain that the constant value is: 



sinh [q(Q = 0)] = - . 

A 



(26) 



From the right panel of Fig. we see that the rotation 
angle 2(p changes by 180° across both the optical-spring 
resonant frequency and the optical resonant frequency. 
The above features in 2ip are typical of resonators and 
can be explained easily from Eq. (|20|) . 

For conventional interferometers (A = 0; continuous 
curves in Fig.^l, the squeeze factor q becomes larger as 
n decreases, providing the strongest squeezing at almost 
all frequencies. In particular, q — > +00 when O — > 0, as 
we can see from Eq. (|26|l . 7 The rotation angle 2(p changes 
by 180° only once over the entire frequency band. 

In the low-power limit (I c — + 0, such that q — ► 0), the 



This gives an upperbound for the amount of squeezing 
achievable with a given optical power, regardless of reso- 
nant features. As we can see from Eq. ()25(l . even for Ad- 
vanced LIGO optical power, at frequencies larger than 
~ 300 Hz, the intrinsic ponderomotive squeezing is al- 
ready very small. 

From the left panel in Fig. Q] we observe that the 
squeeze factor is amplified significantly near the "opti- 



6 In detuned SR interferometers there are two resonances in the 
GW band. One is near the free optical resonant frequency of 
a SR interferometer with fixed mirrors, and we shall denote it 
"optical resonant frequency". The other is shifted up from the 
free pendulum frequency (bel ow 10 Hz) into the detection band 
by the "optical spring" effect l3ll . We shall call it the "optical- 
spring resonant frequency", or the "opto-mechanical resonant 
frequency" . 

7 In reality, q increases only until the test-mass-mirror pendulum 
frequency is reached. 
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transfer matrix M reduces to the rotation matrix 

low-power — R,(7T/2 ^low-power) 5 

with 

2Ae 



Quantity 



Symbol 



/7T 

tan y — - 2<p iow _ pc 



fl 2 - A 2 + e 2 



(27) 



(28) 



Note that this low-power approximation also applies to 
high frequencies where ponderomotive squeezing is sup- 
pressed, even when the power is the typical high power 
in Advanced LIGO [see Fig. [Hand Eqs. |g3J|]. 



B. Rotation of signal quadrature 



Now suppose the output quadrature 
= bi sin C, + hi cos C 



(29) 



is measured, then the signal part of [second term inside 
the parenthesis on the RHS of Eq. is 

S£ oc D\ sinC+D2 cosC oc —A sin£+(e— iVl) cos£. (30) 

Taking the magnitude squared of the above equation, we 
obtain the signal power in this quadrature, 

oc (ft 2 + A 2 + e 2 ) + (ft 2 - A 2 + e 2 ) cos 2( - 2A e sin 2( 

= S + 5iCOs2(C-Cmax), (31) 

where 

5 = (ft 2 + A 2 + e 2 ), (32) 

51 = y/(Q 2 + \ 2 +e 2 ) 2 -4X 2 ft 2 , (33) 



and 



c 



-2A< 



l;i lx --arctan fi2 - - % + ^ 



(34) 



is the quadrature with maximun signal power. Note that 
the relative signal strengths in different quadratures de- 
pend on e and A (i.e., on the optical properties of the in- 
terferometer), but not on i c (i.e., on the laser power and 
mirror masses). Equation (|34() suggests a resonant fea- 
ture of Cmax near the optical resonant frequency, ft ~ |A| . 
Equations H31|) - (|33|l also show clearly a known result [3l|: 
if A ft 7^ 0, we have So > Si, so it is impossible to have 
an output quadrature with no signal. 

By comparing Eqs. (|28|) and l|34|) . we can relate the 
frequency dependence of the maximal-signal quadrature 
of an interferometer with arbitrary optical power to 
the noise-quadrature rotation of the corresponding low- 
power interferometer, that is 



Cn 



1 /7T 

"2 V2 



- 2coi 



(35) 



As we shall see in Sec. IIIII the factor of 1/2 in front of 
the RHS of the above equation makes it difficult to design 
optimal FD schemes near the optical resonant frequency. 



Filter Length 
Input-mirror power transmissivity 

and reflectivity 
End-mirror power transmissivity 

and reflectivity 
Resonant (sideband) frequency 
Bandwidth 



Ti, Ri = 1 — Ti 



7/ 



TABLE II: Parameters of an optical filter. 



C. Noise Spectral Density 

Assuming that ordinary vacuum enters the input port, 
the SR noise spectral density in the C, quadrature (|29|l is 
given by [see e.g., Eqs. (22)-(36) in Ref. 0] 



5'/, 



h 2 

"■SQL 



\Dx sm( + D 2 cos C|' 



[(Cu sinC + C 2 i cosC) 2 + (C12 sinC + C 22 cosC) 2 ] . 

(36) 

By expressing Sh in terms of the quantities </?, Eq. 120( 1. 
q, Eq. 122), Cmax, Eq. ®, and S ,i, Eqs. © and fift, 
we obtain: 



Sh _ , M |2 [cosh 2q - sinh 2q cos 2 (C - y)] ^ 



4tL c n 2 [5 +5lCOS2(C-Crr li «)] 



SQL 



(37) 



III. FREQUENCY DEPENDENT 
INPUT-OUTPUT OPTICS USING KLMTV 
FILTERS 

To realize FD homodyne detection and generate 
squeezed vacuum with FD squeeze angle KLMTV [2(j 
proposed to use Fabry-Perot cavities, detuned from the 
laser frequency, with a transmissive input mirror and a 
perfectly reflective end mirror (ideal case). Later on. 
PC [2l| derived the most general form of the FD quadra- 
ture rotation achievable by these filters. We review their 
work briefly in this section. 



A. Ideal KLMTV filters 

As shown by PC, the most general quadrature rotation 
that can be achieved by a sequence of n ideal KLMTV 
filters, followed by a frequency independent rotation, is 
of the form [see Appendix A of Ref. [2l|]: 



tanC(^) = 



\A n + iB n \>0. (38) 



The complex resonant frequencies of the filters, lu + 
J = 1,2, . ..,71, are given by the roots (with negative 
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fl(£2) 

Input Mirror 



n(Q) 



End Mirror 



FIG. 2: Filter cavity with input a(fl), output b(Q.) and noise 
n(Q) field operators. We indicate with Ti and T e the power 
transmissivity of the input and end mirrors, respectively. 



imaginary parts) of the characteristic equation 

n 

- iB k )Qf = , 3(fij)<0. 



fe=0 



The constant rotation angle is 



9 = arg(A„ + iB n ) . 



(39) 



(40) 



[Our Eq. is different from Eq. (A13) of PC, because 
our definition for flj is different from PC's definition for 
ljj. See their Eq. (A12).] Like in the input-output re- 
lation of SR interferometers, the filter input-output rela- 
tion in this section has also been obtained at the leading 
order in flL/c (as well as in |fi rcs |L/c), that is, in the 
short-filter approximation. It is only under this approxi- 
mation that we can cast the quadrature rotation of these 
filters into the elegant form Ij38(l . 

For low-power interferometers (i c ,q — 0), the trans- 
fer matrix M. reduces to the pure rotation H(ir/2 — 
2<Aow- P owcr) with 



tan 



fix 

^^low-pc 



2Ae 



n 2 



A 2 + e 2 



(41) 



which is of the form 1381) and can be realized by one 
KLMTV filter with complex resonant frequency at luq — 
A — ie, which coincides with the free optical resonant 
frequency of the SR interferometer [see Eq. J7J]. Un- 
fortunately, due to the factor of 1/2 in front of arctan 
in Eq. I|34|) . the frequency-dependent rotation of the 
maximal-signal quadrature cannot be realized by a se- 
quence of KLMTV filters. 



B. KLMTV filters with low loss 

Following KLMTV, we model losses in a filter cavity 
by assuming that the end mirror has a non-vanishing 
power transmissivity T e , and a power reflectivity of R e = 
1 — T e . Denoting the front-mirror power transmissivity 
and reflectivity by T; and Ri (Ti + Ri = 1), the filter 
input-output relation to first order in T e /Ti reads: 



fc out 

/.out 
°2 



= VT^R(C) 



■y/E 



-.filter 

n 

, filter 



(42) 



where the rotation R(£) is the same as in the lossless 
case, nf 1 ^ are vacuum quadrature fields leaking in from 
the end mirror, and the loss factor £ is given by 



'■IE*. 



(43) 



with 



£+ = 



2T e 



i + (±n- W/ )Vr? 



U)f - i*y f 



(44) 

Here ujf, 7/, Lf are the resonant frequency, bandwidth 
and length of the filter, respectively. The bandwidth 7/ 
is related to Tj and Lf by 



7/ 



4L f 



(45) 



[The optical-filter parameters are summarized in Ta- 
ble HH] For a sequence of multiple filters, the rotation 
angles £ and loss factors £ of each filter add up to give 
the total rotation angle and loss factor. In this way, the 
total rotation angle will be identical to the ideal value, 
while the total loss factor will be 



(46) 



The total loss factor is frequency dependent, but never 
exceeds the upper limit 



£n 



E 

J— filters 



42/ 



(47) 



Moreover, if the filters have eigenfrequencies well- 
separated from each other, that is 

\m-\u J f \\ >max{7^7/}, (48) 
and if all filters have high "quality factors," that is 



Luj » 7/ . 



(49) 



then, if we evaluate Eq. (|46|l around the resonant frequen- 
cies only one term dominates, yet away from resonances 
the loss factor is not very large. The total loss factor has 
peaks at the resonant frequencies of each filter, with peak 
value 



23£ 

rpj ' 



j 



(50) 



and width comparable to -fj . 

Once a filter's bandwidth 7/ and the end-mirror trans- 
missivity T e are specified, we can rewrite the peak value 
of the total loss factor (near this filter's resonant fre- 
quency) as 



£ J = 



T 



(51) 
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L f 


7//(27T) 


c 

£-rcs 


4000 m 


100 Hz 


0.0012 




25 Hz 


0.0048 


400 m 


100 Hz 


0.012 




25 Hz 


0.048 


30 m 


100 Hz 


0.16 




25 Hz 


0.64 



TABLE III: Peak values of the filter power-loss factor, £ res 
[Eq. 151H 1. for various filter lengths and bandwidths, assuming 
an end-mirror transmissivity of 20 ppm. 



Thus, the shorter the cavity, the lower the front-mirror 
transmissivity and the larger the loss factor. As an order- 
of-magnitude estimate, we show in Table IIIII the values 
of £ ICB evaluated for typical filter lengths (4000 m, 400 m, 
and 30 m) and bandwidths 7/ (2tt x 100 Hz and 2ir x 
25 Hz) , having assumed T e = 20 ppm [33| . 



C. KLMTV filters with significant loss 

As we can see from Table llHl when the filters are short, 
e.g., on the order of 30 m, the energy loss factor can be- 
come quite large, and the leading-order calculation used 
in Sec. IIIIBl can no longer be trusted. Instead, here we 
give the exact filter input-output relation. By denoting 
with a(O), 6(0) and ro(O) the (Fourier domain) annihi- 
lation operators of the input, output and noise fields at 
frequency iv + Q, we have [see Fig. 0: 



6(0) 



fjT g2i(Q-w / )i / /c 



1 



^T^ e 2i(Q-u: f )L f /c 
^JYjT e i(n- Uf )L f /c 

1 - ^fWRee^-^s^f' 



■a(O) 
n(O) . 



(52) 



Here loq + uif is the resonant frequency of the filter cav- 
ity (the one nearest ujq). The quadrature input-output 
relation can be obtained from Eq. 152|) by using, e.g., 
Eqs. (A8) and (A9) of Ref. 30]. Namely, the relation 



6(±0) = /±(0)a(±0) 



(53) 



valid for annihilation operators is equivalent to the rela- 
tion 



If (/++/*) *(/+-/-) 
2 I -*(/+ - /-) (/+ + /-) 



01 
a 2 



(54) 



valid for quadrature fields 0,1.2 and 61.2- [Note the typo 
in the (2,1) component of Eq! (A9) of Ref. |3Cj .] 

Again, we can apply the short-filter approximation, 
O.L/c, lujL/c, jfL/c <C 1, T e ;$ Ti, and obtain simpler 
formulas: 



1 - a f 



.o 



UJf 



If 



l + Ctf 



.0 



if 



l + a 



.0 



n, (55) 



/ 



7/ 



where a/ = T e /Tj. By converting into the quadrature 
representation we have: 



(^+-/A)(:;)+2^(^) 



l + a f 



.0 



OLf 



UJf 



If 



where ni,2 are quadratures of the field n, 
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l-aj 



2u>/ 
7/ 



If 



2 , .2 



1 - ai H j 



(56) 



(57) 



A 



20 

7 



(58) 



and 



i + a i^L at 

7/ If 



-!£1 

7/ 



If 

1 + a f 



(59) 



IV. SQUEEZED-INPUT AND 
VARIATIONAL-OUTPUT SIGNAL RECYCLED 
INTERFEROMETERS 



A. Input-output relation and noise spectral density 

As discussed by KLMTV, a GW interferometer with 
squeezed vacuum state \S(r,a)) fed into its input port 
can be described by applying the following unitary trans- 
formation. 



«i 

0,2 

\S(r,a) a ) 



R(a)S(r) 

|0a>, 



02 



(60) 
(61) 



in which the quadrature operators undergo a linear 
transformation, while the quantum state is transformed 
back to the vacuum state. [Note that (0a|aj S]-|0s) = 
2-7T Sij <5(0 — 0').] Equation (|6Tj|l suggests that, in practice, 
the squeeze angle a of the squeezed vacuum injected into 
the input port can be obtained by a quadrature-rotating 
optical element, e.g., a KLMTV filter, placed between 
the squeezer and the interferometer. 

Once the unitary transformation is applied, the input- 
output relation of the interferometer can be written sim- 
ilarly to Eq. QJ as 
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FIG. 3: Noise curves of SR interferometers with frequency independent input squeezing and homodyne readout. In the left 
panel, we show the broadband configuration: A = 2n x 234.07 Hz, e = 2ir x 70.36 Hz, £ = —0.8037, with no input squeezing 
(r = 0, a = 0) (continuous curve), with e 2r = 10, a = tv/8 (dashed curve) and with e 2r = 10, a = 37r/4 (dash-dot curve). In 
the right panel, we plot the narrowband configuration: A = 2tt x 600 Hz, e — 2n x 25 Hz, £ = 7r/2, with no input squeezing 
(r = 0, a = 0) (continuous curve), with e 2r = 10, a = 7r/2 (dashed curve) and with e 2r = 10, a = (dash-dot curve). In both 
configurations we fix I c = 840 kW and m — 40 kg and show the SQL curve (dotted line) . 



where 



°11 °12 
°21 °22 



1 

M 



°11 °12 
°21 °22 



e+ r a 2 



Di 

D 2 ) h S QL 



Cix cos a + C12 sin a C\% cos a — Cn sin a 
C21 cos a + C22 sin a C22 cos a — C%\ sin a 



(62) 



(63) 



The quadrature a\ is generally called the "squeezed quadrature" because it enters Eq. I|62(l multiplied by e~ r , while a-i 
is called the "stretched quadrature" because it is multiplied by e +r . If the output quadrature = b\ sin ( + b 2 cos £ 
is measured, the noise spectral density is 



S,, 



er 2r (Cfi sinC + Cft cosCr + ^ (Cf 2 sinC + Cg, cosC)' , 2 



|Z?i sin C + £> 2 cosCI' 



'SQL ' 



(64) 



which in terms of the (ponderomotive) squeeze factor q, intrinsic rotation angle ip, maximal-signal quadrature £ n 
reads: 



S h = \Mf 



[coshq cos(a + C — 293) — sinhg cos(a — ()] 2 + e 2r [coshg sin(a + ( — 2<p) — sinhq sin(a — £)]" 

4et c tt 2 + cos2(C-Cma X )] 

I 



h 



SQL ■ 

M 



In Eqs. and (ffTK)) . the spectral density Sh contains a 
term proportional to e~ 2r , as well as one proportional to 
e 2r . We can take advantage of squeezed vacuum only if 
contains very little (preferably none) of the stretched 
quadrature et 2 . 

In Fig.|3we plot some examples of noise spectral den- 
sities with frequency independent input squeezing (con- 
stant a) and readout (constant £). In this case, squeez- 
ing can improve the sensitiviy at some frequencies, but 
at the price of deteriorating the sensitivity at other fre- 
quencies. However, as investigated by Corbitt, Mavalvala 



and Whitcomb [23L l3~4| , without introducing FD input- 
output techniques, it is still possible to take advantage 
of input squeezing, by choosing carefully the SR parame- 
ters (A, e, l c )i and/or by filtering out the squeezed vacuum 
in the frequency region where the stretched quadrature 
increases the noise. On the other hand, if, for a substan- 
tially detuned configuration, we would like to obtain a 
large noise-suppression factor over the entire frequency 
band, FD input-output techniques should be used. 
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B. Cancellation of the stretched quadrature and 
sub-optimal schemes 

In order that Sh in Eq. (|64|) has only the term propor- 
tional to e~ 2r , we have to impose 

Cf 2 sin C + C 2 Q 2 cos C = , (66) 

or, more symmetrically in a and £, 



f C n 


Cu \ 


( — sin a \ 


{ c 21 


C 22 / 


I cos a J 



It is interesting to note that Eq. (|BTjl does not depend 
on r. This happens because the way e~ r a\ and e r a 2 are 
mapped into [see Eq. (l6*2l ] depends only on a, Cy, 
and but not on r. 



Equation l|67J) can be satisfied in many ways. However, 
since Cy are frequency dependent, either a or £, or both, 
will have to be frequency dependent. Given such a pair 
of (a(Sl), C(^))i the noise spectrum can be obtained by 
inserting them into Eq. I|64|) . obtaining 



-2r 



s h 



C n ( ° 3 sinC(fi) + C 2 T' cosC(^) 



„2r 



C 1 a 2 (f2) sin ((ft) + C 2 a 2 (n) cosC(n) 



h 2 

"■SQL ' 



|£>isinC(fi) + -D2 cosC(^)| 2 

with the second term in the numerator vanishing once Eq. (|67|l is imposed. As a consequence, we can also write 



(68) 



5/, 



-2r 



C n (n) sin ((SI) + cos ((Si) + C^ il > sin C(O) + C 22 w cos ((Si) 



sinC(^) +£> 2 cosC(^)| 2 



-h 2 

"SQL 



_ 2 . r [CiisinC(») + C 2 iC0sC(»)] + [Cu sin((Sl) + C 2 2 cosC(n)] _ 2 
~~ 6 |DisinC(r!)+L'2COsC(fi)| 2 SQL ' 



r 



(69) 



The first equality in Eq. I|69[l says that the noise spectrum 
of an input-output scheme [as specified by (a(Sl),((Sl))] 
with an input squeeze factor r scales as e _2r ; the second 
equality in Eq. I|69[) must hold since for ordinary vacuum 
a rotation of the input quadratures leaves the system in- 
variant. The spectral density, as given by Eq. l|6^|) . is 
e~ 2r times that of a (non-squeezed) FD readout scheme 
with homodyne phase C(^)- Clearly, an additional op- 
timization in ( will give the fully optimal input-output 
scheme. However, we postpone the discussion of the fully 
optimal scheme till Sec. IIVDI and investigate first the 
sub-optimal schemes, which have (a(Sl), C(^)) satisfying 
Eq. (|67|l but do not necessarily have the optimal C(^) 
required by the minimization of (|69|l . These schemes all 
provide a global noise suppression by the factor e~ 2r . 

The (two) simplest solutions to Eq. l(6T|) can be ob- 
tained by imposing ( (or a) to be frequency indepen- 
dent and solving Eq. (|b"7Jl for a (or (). This means that 
KLMTV filters are placed either in the input port or in 
the output port, but not in both places. 

The first simple solution has been studied by Harms 
et al. |25| . who proposed to inject squeezed vacuum with 
FD squeeze angle into SR interferometers. Imposing a 
frequency independent (, they obtained 



tana subop t(^) 



C22 cos ( + C12 sin ( 
C 2 i cos( + Cn sinC 



(70) 



Remarkably, the required a su b pt in Eq. l|7Ufl is of the 
form l|38|) . thus realizable by KLMTV filters. In our no- 



tations, the characteristic equation for the filters is 



Sl 2 (Sl + X - ie)(Sl - X + ie) + [X - 2ie 



ecos 



c] 



= 0, 

(71) 

while the constant rotation following the filters should be 



7T/2-C- 



(72) 



[See Eqs. and (gOJ).] Note that, without mak- 

ing the short-arm and short-filter approximations, both 
Eqs. (|70|l and i)38|) would have been much more compli- 
cated, making the identification of filter parameters much 
less straightforward (or even impossible). 

In this paper, we explore the second simple solution. 
We assume a frequency independent a and requires the 
FD detection phase 



tan ( 



subopt 



(SI) 



C22 cos a — C21 sin a 
C\2 cos a — Cn sin a 



(73) 



This detection phase is also of the form (|38|l and realiz- 
able by KLMTV filters, with characteristic equation 



Sl 2 (SL+X-ie)(Sl-X+ie) + [X + 2e 



esina i r 



0, (74) 



and a subsequent frequency independent rotation 
6 = 3tt/2- a. 



(75) 



Henceforth, we shall call this scheme the BC scheme. The 
noise spectral density of the BC scheme can be obtained 
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by inserting Eq. (|73|) into Eq. I|69|) : the result is Additional insight into these sub-optimal schemes can 

be obtained by decomposing the input-output a-b rela- 
e 2r |-<Wl 2 ^sQL tion into a product of rotation and squeezing operators 

h = 4eL c W\\cosa + (e-in) sina| 2 ' ( ' [see Eqs. ®, © and ©]: 



Afc 



C 



1 ) R(C) 



R(tt/2 - p)S(g)R(-¥>) 



R(a)S(r) 



«i 

«2 



(77) 



readout interferometer input 

Here Ab^ is the fluctuating (noise) part of b^ . Equation l|t)7|) can then be put into the following form: 





1 ) R(C) R(tt/2 - ip) S(q)R(-<p) R(a) 



0. 



(78) 



readout 



interferometer 



input 



r 



In the Harms et al. scheme, the input quadratures 
are rotated (with FD angle a), before entering the in- 
terferometer, in such a way that, after being rotated 
again and ponderomotively squeezed by the interferome- 
ter opto-mechanical dynamics, the squeezed quadrature 
is mapped into a frequency independent output quadra- 
ture, which is detected. In the BC scheme a frequency in- 
dependent squeezed state enters the interferometer. Due 
to rotation and ponderomotive squeezing inside the in- 
terferometer, the squeezed quadrature is mapped into a 
FD output quadrature. We then apply a rotation to the 
field emerging from the interferometer to counteract this 
effect and bring the (image of the) input squeeze quadra- 
ture back to a frequency independent quadrature and 
detect it. 

Finally, another interesting sub-optimal scheme can be 
obtained by imposing £ = a = —tp. In this case the noise 
part of the output quadrature field l|77l) is 



(0 l)R(7r/2)S(g)S(r) 



a i 

0,2 



e-^+^a! (79) 



which gives the lowest amount of noise (but does not 
guarantee a maximal signal content). Unfortunately, 
from Eq. I|21|) we see that tan ( = — tan ip is not of the 
form (J3SJ, and thus not realizable by KLMTV filters. 



C. Sub-optimal schemes using q— ip 
parametrization: the low-power limit 

If the ponderomotive squeezing factor q is small, the 
fully optimal input-output scheme can be solved easily 
using the various quadrature-rotation angles. As seen 
in Sec. [HJ a small q can either arise from a low optical 
power, or from considering high frequencies (/ si 300 Hz 
for Advanced LIGO power), see Eq. (|2ip. However, we 



shall still refer to this as the low-power limit. In this 
case, the output noise is proportional to 



a i 

«2 



(0 1)R(C + tt/2-2 

^low-power T" OL) ) 

and the minimal noise is obtained whenever 

Q ^ Ct 2 ^low-power ■ 

By setting £ equal to the maximal-signal quadrature [see 
Eq. (SEN, 



(80) 



(81) 



Cnaax " ^ 



ow- power j 



(82) 



we find the fully optimal readout scheme: 



(Ca)opt = (-4 



^low-power j A ~r V^low-powcr I ■ (83) 



Simple as it looks, this fully optimal scheme is not re- 
alizable by KLMTV filters because tan £ opt and tan a opt 
given by Eq. (H21) are not of the form (fTT^fl . 

We now compare the Harms et al. (H) and BC schemes 
in the small-g regime. They can be written in terms of 
(£, a) as 



(C) a )n — (C) 2 ^low-power — C) > 
(C) a )BC = (2 ^low-power Ot, Ct) ■ 



(84) 



The two schemes give the same noise output part e r a\ 
[see Eq. (I80(l ], while for the signal power they yield [see 
Eqs. (EUl-jSSl)] 

S H = S + Si COS [2(C - Cmax)] 

= S + St C0S(2 C + 7t/2 - 2 ^low-power) , (85) 



and 



= So + Si COS [2(C - Cmax)] 

= So + Si COs(-2 a + 7r/2 + 2 <£>low-power) 

= So + Si cos(2 a - 7r/2 - 2 ^i w- P owor) ■ (86) 
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FIG. 4: Equivalence of the Harms et al. and BC schemes at high frequencies. We assume J c = 840 kW, m = 40 kg and e 2r = 10. 
In the left panel, we plot the broadband configurations with A = 2n x 234.07 Hz, e — 2n x 70.36 Hz, while in the right panel we 
show narrowband configurations with A = 2ir x 600 Hz and e — 2ir x 25 Hz. In both panels, BC schemes with a = ir and 7r/2 are 
shown in dark continuous and dark dashed curves, respectively, while Harms et al. schemes with f = tt/2 and £ = are shown 
in light continuous and light dashed curves. Noise curves from the two schemes, which are related by a — £ + 7r/2, do agree 
quite well for frequencies higher than 200 Hz in the broadband configuration (left panel) , and at almost all frequencies in the 
narrowband configuration (right panel). Note that at high frequencies q is significantly lower in the narrowband configuration, 
see left panel of Fig. [3] 



This means, the two sub-optimal schemes have the same 
ideal performance in the low-power regime and we can 
map one into the other by setting a <-> £ + 7r/2. 

This equivalence can be understood more intuitively 
if we compare the dependence of the various readout 
quadratures (i.e., maximum-signal, Harms et al., and 

BC) On ^low-power- 

The maximal-signal quadrature Cmax rotates as const + 
<y5iow-powcr ■ In the Harms et al. scheme, the detected 
quadrature is constant, and therefore lags the maximal- 
signal quadrature by const + yiow-power- In the BC 
scheme, the detected quadrature rotates as ( = const + 
2<Aow-powcr, which advances the optimal quadrature by 
const + yiow-power- In this way, if one adjusts the con- 
stants (by adjusting £ in the Harms et al. scheme and a 
in the BC scheme), the detected quadratures in the two 
schemes can be made to lie symmetrically on each side of 
the maximal-signal quadrature. Since the detected sig- 
nal power depends only on cos [2(£ — (max)] [see Eq. IpTTjl]. 
which is an even function of (C — (max), the two schemes 
must detect the same signal power and hence have the 
same sensitivity. 

In Fig.QJ we give examples of the BC and Harms et al. 
noise curves for two SR interferometers, a broadband con- 
figuration (with A = 2?r x 234.07 Hz, e = 2tt x 70.36 Hz) 
and a narrowband one (with A = 2ir x 600 Hz and 
e = 2ir x 25 Hz). For both interferometers, we use 
I c = 840 kW, m — 40 kg. Although the optical power 
used here is not low by any practical standards, the two 
schemes for the broadband configuration already agree 



quite well, under the correspondence a <-> £ + tt/2, for 
frequencies above ~200Hz. The two schemes are equiva- 
lent for the narrowband configuration for almost all fre- 
quencies. The better agreement in the narrowband con- 
figuration can be understood easily by realizing that pon- 
deromotive squeezing is weaker in this case, as shown in 
the left panel of Fig. [3] 

D. The fully optimal scheme and the BC scheme at 
low frequencies 

In this section, we consider the fully optimal scheme. 
Analytical formulas of the fully optimal detection 
quadrature has been obtained by Harms et al., but we 
provide an alternative approach, yielding results in sim- 
pler form and more related to the BC scheme. 

It is straightforward to show that (as also done by 
Harms et al. and reviewed in Appendix lA"|l . fixing £ and 
r, the a obtained from Eq. I|70|l gives the (constrained) 
minimum noise. On the contrary, fixing a and r, the 
readout quadrature C obtained from Eq. I|73l) does not 
give the constrained minimum. Instead, minimizing Sh 
[Eq. I|64l) ] over £ (with a fixed) requires a rather compli- 
cated readout phase, determined by one of the two roots 
of 

T 2 tan 2 C + T\ tan C + Fq = , (87) 

where 
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FIG. 5: Approaching the fully optimal input-output scheme by taking the envelope of Harms et al. and BC noise curves. We 
fix I c = 840 kW and m = 40 kg, and A = 2vr x 234.07 Hz, e = 2vr x 70.36 Hz, and assume e 2r = 10. In the left panel, we 
plot the noise curves of Harms et al. schemes, with frequency independent readout phase ( = (light short dash), 7r/4 (light 
long dash), n/2 (light dot dash) and 37r/4 (light dot), and FD input squeeze angle given by Eq. 17011 . These curves cross each 
other near both the optical resonance and the optical-spring resonance. The fully optimal curve (dark continuous) is obtained 
by taking the lower envelope of the entire family of these curves. In the right panel, we plot the BC curves, with frequency 
independent input squeeze angle a = — n/2 (light short dash), — n/4 (light long dash), (light dot dash) and n/4 (light dot), 
and FD readout phase given by Eq. H73II . The lower envelope of these curves also gives the fully optimal noise curve (dark 
continuous, indentical to the one shown in left panel). The fact that these curves do not cross each other at low frequencies 
suggests that one member of this family is fully optimal in this band. Indeed, the BC curve optimized for Q — [dark dashed 
curve, with a = arctane/A, see Eq. 1)91^ ] does agree very well with the fully optimal curve for frequencies lower than ~ 200 Hz. 
In both panels we also show the SQL line. 



?2 = [(C?r r ) 2 + (C?f) 2 } K(DID 2 ) - (C?{- r C 2 r r + CVT C£ r ) l^l 2 , (88a) 

*i = ((c?0 2 + (cvY) 2 ) \Di\ 2 - mcf + (c 2 Y) 2 ] |A| 2 , (88b) 

n = (err c 2 °r r + cs r c 2 t) \d 2 \ 2 ~ [(c«o 2 + (c 2 Y) 2 ] r^m • 



Equations (|87|l - H88c[l . which we obtained indepen- 
dently |35| from Harms et al., are equivalent to Eqs. (28)- 
(30) of Harms et al. once we set r to zero in Eqs. Q88a[l - 

As said above, the fully optimal scheme, denoted by 
(a op t(fi), Copt(^))) should satisfy the sub-optimal condi- 
tion (|67|l . As a consequence, the noise spectrum of the 
fully optimal scheme is also given by Eq. (|69|l . when C(^) 
is replaced by £ opt (f2). Therefore, Copt(^) can be ob- 
tained by minimizing the Sh in Eq. (|69|) . which is given 
by the special case of Eqs. l|S7j) - (|88c|) with r = [or 
Eqs. (28)-(30) of Ref. a opt (ti) can then be obtained 
from Eq. (J7DJ). It is evident from Eq. JS7JI that the fully 
optimal scheme cannot be realized by KLMTV filters, 
except in special cases, e.g., for conventional interferom- 
eters. As observed by Harms et al., the optimal noise 
spectrum can also be obtained graphically, by plotting 
all the noise curves with different constant values of £, 



and then taking the lower envelope of all these curves, as 
seen in the left panel of Fig. (and Fig. 4 of Ref. |2jJ). 
The optimal £ at each frequency is the one whose noise 
curve touches the envelope. 

We now deduce the optimal scheme in another way. 
Again, since (a op t(0), Copt(^)) satisfy Eq. l|rJ7ll . the fully 
optimal noise spectral density can also be obtained by 
taking the minimum among all BC noise spectral den- 
sities with all possible a — the minimum is achieved 
automatically in a op t(f2), and for Copt ( ^ ) it is given by 
Eq. I|73(l . Similarly, this can be done graphically by tak- 
ing the lower envelope of BC noise curves with all pos- 
sible a, as shown in the right panel of Fig. |3J From 
the plot, it is interesting to observe that, there are no 
crossings between different BC noise curves at low fre- 
quencies (differently from the Harms et al. curves in the 
left panel) , suggesting that one BC curve might be nearly 
fully optimal at these frequencies! 
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More quantitatively, since the BC noise spectrum l(7B|) see Eqs. (|8"Tj) - H88cl) ]: 
has a much simpler dependence on a (than the depen- 
dence of the Harms et al. noise spectrum on Q, it is t 2 ^ £ (89) 
much simpler to obtain the optimal input squeeze angle opt — X 2 _ e 2 _ f]2 ' v > 
a op t from this approach [than to obtain £ opt from the 
approach starting with the Harms et al. noise spectrum, and 



gopt 



\M\ 2 hl 



SQL 



2 et c ft (A 2 + e 2 + Q 2 ) 



(90) 



2\n 



A 2 



Q 2 



These simple explicit expressions of >!5 opt and a op t(ft) 
have not been previously obtained. The optimal read- 
out phase Copt can be obtained from Eq. (|73[) . From 
Eq. (|89|l . we can see that the fully optimal scheme can- 
not be achieved by KLMTV filters. The only exception 
is when A = (i.e., for a conventional interferometer). In 
this case we have a op t = 0, and the £ pt(ft) is given by 
Eq. and it is reliazable by KLMTV filters. This is 
exactly the KLMTV squeezed- variational scheme. 

Although the form of a pt is not achievable by KLMTV 
filters, we note that, at low frequencies (lower than the 
optical resonant frequency), the variation in a opt is mild. 
In fact, by setting in the BC scheme 



ot = a op t(ft = 0) = arctan \ — j , 



(91) 



we obtain . 



oBC low-frcq 



e- 2r \M\ 2 h 2 



SQL 



2 e i c n 2 



2(A 2 + e 2 + fl 2 ) - 



2\ 2 n 2 



x 2 



(92) 



Taking the ratio between S? C low frcq and «S^ pt , and ex- 
panding in f2, we have 

-rBC low-frcq / \ \ 2 

= 1+(J^] f— ! . (93) 



gopt 



A 2 + e 2 



A 2 



The correction factor in Eq. I|93|l is usually small at low 
frequencies. For example, by maximizing over either e or 
A, it is easy to show that 



Ac 



n 2 



A 2 



< 



27 



(94) 



256 \max{A, e} / 

at worst. The correction in the noise spectral density can- 
not exceed ~ 10% (in power) for Q ~ max{A, e}. For sub- 
stantially detuned configurations (A exceeding ~ 200 Hz), 
this makes the BC scheme essentially fully optimal up to 
~ 200 Hz. This result is confirmed by the right panel 
of Fig. [31 in which S^ C low flcq is plotted (dark dashed 
curve) in comparison with (dark continuous curve). 



V. APPLICATIONS TO ADVANCED LIGO 

In this section, we discuss the possibility of applying 
the above FD techniques to Advanced LIGO interfer- 
ometers. As shown by KLMTV 20], a major difficulty 
in making those techniques practical for advanced inter- 
ferometers is the issue of optical losses. Given a cer- 
tain bandwidth and mirror quality (i.e., round-trip loss 
in the filter cavities), the shorter the filters, the higher 
their optical losses (see Table ITlT|) . In fact, in order to 
achieve third-generation performance, optical filters in 
the squeezed-variational scheme will have to be ~ kilo- 
meter in lengths. In Advanced LIGO, kilometer-scale fil- 
ter cavities are not practical and only short filters can fit 
into the corner-station building. A plausible length scale 
is ~ 30 meters; and the realistic round-trip loss is around 
20ppm [33. With such short (and lossy) filters, we shall 
assume most of the time that filter losses will dominate 
and ignore internal interferometer losses [see Sec. V in 
Ref. [31| for treatment of lossy SR interferometers] . We 
shall only comment briefly on the effect of internal losses 
when discussing narrowband sources. The noise spec- 
trum with filter losses are obtained by using the exact 
input-outp ut rel a tion o f KL MTV filters (Sec. HUE). 

In Sees. IV Al IV Bl and IV CI respectively, we shall 
discuss the broadband configuration optimized for the 
detection of NS-NS binary inspiral waveforms, the nar- 
rowband configuration targeting GWs from specific ac- 
creting NS's and the wideband configuration that can be 
used to observe several kind of sources. [For an exhaus- 
tive discussion and summary of GW sources for advanced 
interferometers see, e.g., Ref. [36}.] 



A. Broadband configuration: NS-NS binary 
inspiral 

Inspiral waves from compact binaries (NS-NS, NS- 
BH or BH-BH) are among the most promising sources 
for Advanced LIGO. In this section, we discuss the so- 
called broadband configuration obtained by maximizing 
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Interferometer Configuration 


Filter I 


Filter II 


Performance 


Input-Output Mirror „— 2r e A ^ 
Scheme Type 2tt Hz 2it Hz 


2tt Hz (ppm) 


27rHz (ppm) 


SNR Event Rate 
300 Mpc Improvement 


No Squeezing Spherical 1 70.4 234.1 -0.804 
No Filters 0.1 561.8 55.2 1.522 -0.040 
Harms et al. 0.1 280.4 296.3 FD -0.381 
BC 0.1 157.4 355.2 -2.090 FD 


296.0 - 285. 4t 717 
352.0 - 161. 0i 404 


-59.8 - 23. 7i 59.6 
-59.5- 12. 3i 30.9 


5.44 1.00 
6.73 1.89 
7.15 (7.81) 2.27 (2.96) 
6.77 (7.84) 1.93 (2.99) 


No Squeezing M ^ a " 1 14.4 179.1 -1.010 

No Filters 0.1 275.3 101.3 1.395 -0.131 
Harms et al. 0.1 106.5 233.9 FD -0.518 
BC 0.1 55.7 240.5 -2.179 FD 


227.2 - 115. 4t 290 
230.8 - 61.2 i 154 


-73.2- 18. Oi 45.1 
-72.7- 9. Si 24.6 


9.29 1.00 

10.44 1.42 
11.73 (15.08) 2.01 (4.27) 
10.45 (15.63) 1.42 (4.76) 



TABLE IV: Optimizations of SR interferometers with (i) no squeezing, (ii) frequency independent squeezing and homodyne 
detection ("no filters") (iii) FD squeezing but frequency independent readout (the Harms et al. scheme), and (iv) frequency 
independent squeezing but FD readout (the BC scheme) for neutron-star-binary inspirals, with quantum noise, seismic noise 
and thermoelastic noise (with spherical and Mexican-Hat mirrors) included. The only optical losses included are those from 
the 30-m optical filters. The round-trip loss of each filter is set to be 20 ppm |3^l . Noise curves of configurations listed here are 
plotted in Fig. HJ For the Harms et al. and BC schemes, we also optimize the SNR when there is no optical losses, those SNRs 
and the corresponding event-rate improvements are quoted inside brackets. 




FIG. 6: Noise curves optimized for NS-NS binaries, for no squeezing ("no squeezing", light dashed curves), frequency- 
independent squeezing and homodyne detection ( "no filters" , light continuous curves) , the Harms et al. scheme (dark continuous 
curves) and the BC scheme (dark dashed curves). Parameters of each configuration are listed in Table llVl Quantum noise, 
seismic noise and thermoelastic noise of sapphire (also shown, in dash-dot curves) are included to give the total noise curves. 
We have used predictions for the thermoelastic noise of spherical mirrors (left panel) and that of Mexican-Hat mirrors |3?| 
(right panel). In addition, we have shown the characteristic strengths of possible GWs from LMXBs (diamond, solid circles, 
solid squares and open triangles) and known radio pulsars (thin dashed lines). 



the signal-to-noise ratio for NS-NS inspiral waveforms, 
proportional to 




Sh(f) 



df, 



where 



\kf)\ = ^r 7/6 e(/i 



ISCO 



f) 



(95) 



(96) 



is the frequency-domain amplitude of the leading (New- 
tonian) order inspiral signal in the stationary-phase ap- 
proximation. The cutoff frequency is chosen to be /iscOj 
the GW frequency corresponding to the Innermost Sta- 
ble Circular Orbit (ISCO) of a Schwarzchild black hole 



with mass 2Mns = 2 x 1.4 M = 2.8 M Q , which is equal 
to 1570 Hz. In the optimization, we have also included 
the seismic noise 



S% eis = 5.3 x 10 



-21 



10 Hz 



9.2 



(97) 



and the thermoelastic noise of sapphire mirrors with 
spherical surfaces, as in the baseline design, 



S£ E - 2.7 x 10- 24 



100 Hz 



(98) 
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as well as when the so-called Mexican-Hat mirrors are 
used, which are designed to reduce this noise |37| . 

V^™ = 1.1 x ID- 4 ( -i= . (99) 

In Table IIVI we list the values of A, e, a (frequency 
independent squeezing angle for BC scheme), £ (fre- 
quency independent detection quadrature for Harms et 
al.'s scheme), obtained by optimizing the SNR of NS- 
NS binary inspirals at 300 Mpc, and the corresponding 
optimal SNR. We assume I c = 840 kW, m = 40 kg, 
and e~ 2r = 10 and did the optimization for (i) non- 
squeezed SR interferometers, (ii) SR interferometers with 
frequency independent squeezing and homodyne detec- 
tion ("no filters"), 8 (hi) squeezed SR interferometers with 
the Harms et al. scheme (FD input squeezing + ordinary 
homodyne detection) and (iv) squeezed SR interferome- 
ters with the BC scheme (ordinary squeezing + FD ho- 
modyne detection). In the Table we also give the im- 
provements in the predicted event rate with respect to 
non-squeezed configurations, as the cube of the improve- 
ments in SNR at a fixed distance. 

As we can read from Table IIVI with frequency inde- 
pendent squeezing (i.e., no filters), it is already possible 
to improve the NS-NS event rate by a significant amount, 
89% (spherical mirror) or 42% (MH mirror). The Harms 
et al. scheme provides further improvement in the event 
rate with respect to the no-filter case, by 20% (spherical 
mirror) or 42% (MH mirror). The BC scheme, however, 
being more susceptible to filter optical losses, does not 
yield as good a performance. In order to appreciate how 
much the filter optical losses affect the sensitivity, we 
have also optimized the SNR for the FD schemes without 
including filter optical losses (but with thermal and seis- 
mic noises included), the results are quoted in brackets 
in Table ITVl FD schemes without losses can outperform 
frequency independent squeezing significantly. For exam- 
ple, the ideal BC scheme can have 58% (spherical mirror) 
or 235% (MH mirror) more event rates than the no-filter 
case. [In this case the the BC scheme can also provide 
slightly higher event rates than the Harms et al. scheme, 
due to better sensitivity at low frequencies (but mostly 
still masked by the thermal noise), by 1% (spherical mir- 
ror) or 11% (MH mirror).] 

Noise curves corresponding to the optical configura- 
tions listed in Table IIVI are plotted in Fig. We notice 
that due to optical losses the BC noise spectral densities 
have a peak around the optical-spring resonant frequency. 
The noise spectrum of the "no filters" scheme (squeezing 
with frequency-independent input-output optics) is com- 
parable to the Harms et al. and BC schemes at high fre- 
quencies, but becomes worse at low frequencies. These 
"no-filter" curves are quite similar to the wideband noise 



Corbitt, Mavalvala and Whitcomb l34l are currently investigat- 
ing this scheme. 



curves proposed by Corbitt and Mavalvala [2j| , especially 
in the case of spherical mirrors. 

The squeezing noise curves optimized for NS-NS bi- 
naries also have better high-frequency sensitivity than 
non-squeezed configurations, although they were not op- 
timized specifically for high frequencies. From Fig.® we 
see that for frequencies higher than ^500 Hz (spherical 
mirrors) or ~300 Hz (Mexican-Hat mirrors) , the squeezed 
configurations are ^5-8 (spherical mirrors) or ~ 3 - 5 
(Mexican-Hat mirrors) as sensitive (in amplitude) as the 
non-squeezed configurations. [The Mexican-Hat mirrors 
produce lower thermoelastic noise, so the noise spectral 
densities are better optimized at low frequencies, reduc- 
ing the bandwidth. This is why in this case the noise 
curves optimized for NS-NS binaries yield worse high- 
frequency sensitivity than those with spherical mirrors.] 

In Fig. we also plot (in light thin dashed lines) the 
characteristic GW strengths from known radio pulsars. 
Following the notation of Cutler and Thorne |3y|, the 
characteristic strength St c is defined as the maximum 
allowed noise spectral density Sh(f source) (at and near 
the source frequency / SOU rcc) such that the source is de- 
tectable. Note that Sh c will in general depend on the 
data analysis technique and statistical criteria used, e.g., 
integration time, confidence level, etc.; sometimes it is 
also obtained by averaging over unknown source param- 
eters, such as the spin orientation of pulsars [see App. [5] 
for more details]. Here for known radio pulsars at lOkpc 
distance, with ellipticity e = 10~ 6 , 10~ 7 and 10~ 8 , we 
have been assuming 1% of false-alarm probability in a 
coherent search of 10 7 s of data (coherent search for such 
a long time can only be done for pulsars whose sky posi- 
tions and phase evolutions are known [38, 39]). 

From Fig. [5] we see that the NS-NS optimized noise 
spectra for spherical mirrors can detect known pulsars 
at lOkpc with e > 10~ 7 if the GW frequency is higher 
than 500 Hz, while those for MH mirrors can detect e > 
2 x 10~ 7 if GW frequency is higher than 300 Hz. 

We have also shown in Fig. the frequencies and the 
estimated characteristic GW strengths from LMXBs (Sco 
X-l in diamond, the Z sources in solid dots, Type-I 
bursters in solid squares, and accreting millisecond pul- 
sars in open triangles). We shall explain those sources in 
more detail in the next section and in App. [Bj All the 
squeezed-input configurations are able to detect Sco X-l 
with large margins, while configurations with spherical 
mirrors might also be able to detect the group of six Z 
sources near 600 Hz. 



B. Narrowband configuration: LMXB 

Low-Mass X-ray Binaries (LMXBs) are systems 
formed by a neutron star and a low-mass stellar com- 
panion, from which the neutron star accrets material. 
Observations of LMXBs have provided evidence of a NS 
spin- frequency "locking" in the range 260 Hz < f s < 
600 Hz (much lower than the breaking frequency of ~ 
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0.00 
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0.01 
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0.03 
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370 


0.25 


0.08 


0.01 


0.03 


0.03 
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TABLE V: Frequency, characteristic amplitude h c , and characteristic strengths (y S h ay , 20-day coherent integration with 

1% false alarm) of possible GWs from several known LMXBs [inclu ding Z sources (the first 8 sources), Type-I bursters (the 
next 8 sources) and accreting millisecond pulsars (the last 3 sources)] |44jl. and sensitivities achievable by non-squeezed, Harms 
et al. and BC schemes. Both 5dB (e _2r = 0.316) and 10 dB (e _2r = 0.1) squeezing are considered. Sensitivity is measured 
by taking the ratio between the characteristic strength and square root of the noise spectral density at the predicted GW 
frequency. Bold face is used for ratios larger than unity, in which case the GW is detectable. In this Table we use the baseline 
assumption (al) (i.e., f s = fd f° r Z sources, and mass-quadrupole emission). To obtain the predictions for other mechanisms 
of GW emission and determination of the spin frequency, see Table I VII Note that Sh c should convert by the same factor as 
h c . For the source SAX J1808.4— 3658 in particular, we have also shown in parenthesis values that correspond to a 4-month 
coherent integration. 



1.5 kHz (43). These systems are rather old and believed 
to have been spun up by accretion torque. Thus, to 
explain the locking it has been conjectured that accre- 
tion torque could be balanced by angular-momentum loss 
due to GW emission [U ED, E3- In Table we list 
a number of LMXBs that are promising GW sources: 
the first group contains the so-called Z sources, the sec- 
ond group the Type-I bursters, and the third group 
accreting millisecond pulsars (all data are taken from 

Refs. ei mm). 

The spin frequency of the NS in these LMXBs is not 
unambiguously determined, except for accreting millisec- 
ond pulsars, whose X-ray fluxes pulsate at their spin fre- 
quencies, i.e. fp = f s . For Type-I bursters, the spin fre- 
quency can be inferred from the millisecond oscillations 
in their X-ray fluxes observed after bursts (fs) and from 
the kHz QPO difference frequency (fd)- However, for dif- 
ferent sources, it has been observed that either fd = fp 
or fd = /b/2, and it is not firm yet whether f s should be 
equal to fp or fd- Recently, X-ray bursts have been ob- 
served from the source SAX J1808. 4-3658 (an accret- 
ing millisecond pulsar, with spin frequency known from 
f s = fp |4(| ) , and X-ray flux after the bursts is observed 
to oscillate at the spin frequency (i.e., fp = fp = f s ). 
Moreover, for this source the kHz QPO difference fre- 
quency is observed to be half this value: fd = fp/2. 



1 (MQ) 

2 (CQ) 



a (fs = fd) 



b (fs = 2fd) 



(al) ,(al) 
GW ' " c 



n(«i) 



2 f (al) 
3JGW ' 



(al) 



3 J GW > V 4 c 



TABLE VI: Conversion of predicted GW frequencies and 
characteristic GW amplitudes (/gw,/i c ) from LMXBs be- 
tween different assumptions on spin frequency and GW emis- 
sion mechanism. In particular, our baseline assumption, (al) 
for Z sources and (1) for Type-I bursters and accreting mil- 
lisecond pulsars, has been used by Refs. pSfil I43L El ] to give 
numerical estimates for GW frequency and characteristic am- 
plitudes/strengths. [For Type-I bursters and accreting mil- 
lisecond pulsars, the conversion from (1) to (2) follows the 
rule from (al) to (a2) in the table.] 



This might favor the argument that /, s = fp for all Type- 
I bursters, as assumed by Refs. j^J Hj, |4J] and used in 
Table [V] [henceforth we shall always adopt this assump- 
tion]. For Z sources, only kHz QPOs have been observed; 
this makes it difficult to determine the NS spin frequency: 
it could be either (a) f s = fd or (b) f s = 2fd [note that 
for different Type-I bursters either (a) or (b) could be 
true] . 

Moreover, two plausible physical mechanisms for GW 
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Interferometer Configuration 


Filter I 


Filter II (unapplied) 


Performance 


_ 2r e A 

2irHz 2ttHz 


27T Hz (ppm) 


2ir Hz (ppm) 


V^7(10- 24 /\/H5) BW NS-NS 300 Mpc 
at 600 Hz (Hz) Sphcrical/MH 


No Squeezing 1 (OdB) 25 601.4 -0.748 






0.89 41 2.93/4.43 


Harms et al. 0.316 (5dB) 60 601.2 FD -0.749 
BC 60 601.3 0.806 FD 


599.7 - 60.3 i 152 

599.8 - 60.3 i 152 


-43.0-2.2i 5.6 
-43.0- 2. 3i 5.8 


0.94 152 4.59/6.70 
0.95 153 4.59/6.64 


Harms et al. 0.1 (lOdB) 100 597.9 FD -0.722 
BC 100 598.3 0.843 FD 


596.4 - 100.5 i 253 
596.7 - 100. 5 i 253 


-44.1- 3.8 i 9.7 
-44.1- 3. 9i 9.7 


0.88 356 6.03/8.15 
0.89 361 5.98/7.96 



TABLE VII: Optimization of SR interferometers with (i) no squeezing, (ii) FD squeezing but frequency-independent readout 
(the Harms et al. scheme), and (iii) frequency-independent squeezing but FD readout (the BC scheme) for narrowband sources 
around 600 Hz. We have considered both 5dB (e _2r = 0.316) and 10 dB (e" 2r = 0.1) squeezing. In both of the FD schemes, 
filter II, which has impractically high finesse, does not affect high frequency performance, and is not applied. The 600 Hz 
sensitivity, bandwidth, and SNR for NS-NS binaries at 300 Mpc are given as performance indices. Here bandwidth is defined as 
the difference in the two frequencies at which sj Sh (/) = \/2Sh(600 Hz). Noise curves of configurations listed here are plotted 
in Fig.|71 



emission from accreting NS's have been proposed: (1) 
mass quadrupole radiation from deformed NS crusts 
(/cw = 2/ s ) 0)0]; and (2) current quadrupole radi- 
ation from unstable (with respect to gravitational radi- 
ation) pulsation modes (r-modes) in NS cores (/gw = 
4/ 5 /3) |H E3 EH- Suppose one of the two emis- 
sion mechanisms to dominate, then along with uncer- 
tainties in spin frequencies, we have four possibilities 
for Z sources, (al), (a2), (bl) and (b2); [and two pos- 
sibilities for accreting millisecond pulsars and Type-I 
bursters, (1) and (2)]. In the following, we consider (al) 
for Z sources and (1) for accreting millisecond pulsars 
and Type-I bursters our baseline assumption, as done in 
Refs. [3a, El EH> [!-!■ an d comment on what happens if 
the other options turn out to be true. 

In the second column of Table E] we list GW frequen- 
cies obtained from the baseline assumption; GW frequen- 
cies based on other assumptions can be obtained from 
the (al) or (1) value by using Table IVII The charac- 
teristic GW amplitude h c from LMXBs has been esti- 
mated |48l |49| by assuming a balance between GW an- 
gular momentum loss and accretion torque, with the lat- 
ter estimated from X-ray flux, and by subsequent av- 
eraging over the (unknown) spin orientation [see Ap- 
pendix IB II for a detailed explanation on the averaging 
process and the associated uncertainties]. However, the 
value of h c can also be different due to the various as- 
sumptions on spin frequency and GW emission mecha- 
nism we can make. Values listed in the third column of 
Table|V]has been obtained in Refs. EJESE3 using the 
baseline assumption; conversions from (al) or (1) to the 
other assumptions can be made easily using Eq. (8) of 
Ref. and Eqs. (4.4)-(4.6) of Ref. 49], and are given 
in Table IVll By assuming 1% false-alarm probability and 
20-day coherent integration time [due to unknown orbital 
motion and frequency drifts caused by fluctuations in the 
mass accretion rate] Sh c can be obtained from h c (listed 
on the fourth column of TablelVl see Add. IB 2l for details), 
note that for the different assumptions -J 'Sh e changes by 
the same factor as h c . For the accreting millisecond pul- 
sar SAX J1808. 4-3658, for which the orbital motion is 
known |46j |. assuming that GW frequency evolution can 
be obtained, we also show (in parenthesis) the character- 



assumption 



A 



Hz 2ttHz 2ttHz 



C 



(al) 600 30 600.4 -0.756 

(bl) 1200 90 1057.3 -0.065 

(a2) 400 25 412.6 -0.749 

(b2) 800 90 769.0 -0.306 

TABLE VIII: Parameters of narrowband configurations tuned 
to LMXB sources when different assumptions on spin fre- 
quency and GW emission mechanism are adopted. Noise 
curves of these configurations, with and without interferome- 
ter losses included, are shown in Fig. and compared to the 
corresponding GW characteristic strengths. 



istic strength obtained with a 4-month integration. 

It is important to realize that there are still uncertain- 
ties as to whether a particular source will be detectable, 
even if the noise curve is below Sh c — as explained in 
Appendix [BJ However, the main aim of this paper is 
to discuss interferometer configurations, rather than the 
data analysis of narrowband sources, so we shall use Sh c , 
as done by Cutler and Thorne [3g despite the subtleties, 
as a playground to compare sensitivities of different inter- 
ferometer/filter configurations. Conclusions drawn in our 
discussions on whether these sources will be detectable 
should definitely be refined by more rigorous investiga- 
tions. 

In Fig. we plot the noise curves obtained for a non- 
squeezed SR interferometer and for squeezed SR interfer- 
ometers with the Harms et al. and BC schemes by opti- 
mizing their sensitivities in a narrow band around 600 Hz. 
Peak sensitivities and bandwidths are adjusted to incor- 
porate the signal strengths of a group of 7 Z sources (in- 
cluding ScoX-1). The baseline assumption is used in 
obtaining /gw and Sh c for these sources. 

For the non-squeezed interferometer, we obtain a noise 
curve similar to the "narrowband" curve in Fig. 1 of Cut- 
ler and Thorne |36j |. provided originally by Ken Strain. 
For squeezed interferometers, we have considered both 
5dB (e- 2r = 0.316) and 10 dB (e' 2r = 0.1) squeezing. 
Since in narrowband configurations, the seismic and ther- 
mal noises do not affect significantly the choice of the SR 
parameters, the noise curves in Fig. \7\ have been opti- 
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FIG. 7: Noise curves of non-squeezed (light dashed curve), 
Harms et al. (dark continuous curves) and BC (dark dashed 
curves) configurations optimized for narrowband sources, for 
5dB and 10 dB squeezing. We apply only one filter, that is 
the one with resonant frequency near the free optical resonant 
frequency of the SR interferometer, or filter I (see Table Wl . 
The interferometer noise curves contain only quantum noise 
but include filter losses. The thermoelastic noise of spheri- 
cal and MH mirrors, and the SQL are plotted for comparison. 
We also show the frequencies and characteristic strengths (20- 
day coherent integration, 1% false alarm [3s3,|3tJ) of possible 
GWs from LMXBs [4l| [using the baseline assumption (al), 
namely / s = fd and mass-quadrupole emission]: Z sources in 
solid circles (Sco X-l in diamond), Type-I X-ray bursters in 
solid squares and accreting millisecond pulsars in open trian- 
gles. For the accreting millisecond pulsar SAX J1808. 4-3658, 
for which the orbital parameters and GW phase evolutions 
are known |45l| . we show in another open triangle (linked to 
the 20-day one with a vertical segment of solid line) the char- 
acteristic strength assuming a 4-month integration. 
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FIG. 8: Relative increase in SNR for LMXB sources around 
600 Hz, with 5dB (e" 2r = 0.316, dashed curve) and 10 dB 
(e~ 2r = 0.1, continuous curve) squeezing. Since the Harms et 
al. and BC schemes are extremely close to each other only one 
curve is shown for each squeeze factor. The various detectable 
LMXBs [under the baseline assumption] listed in Table IV! are 
also shown (ScoX— 1 in solid diamond, the rest in solid circles) 



mized using only the quantum-optical noise (but we in- 
clude filter optical losses). [For comparison we plot in 
Fig. the thermoelastic noises.] We obtain the parame- 
ters e, A and £ for the squeezed configurations following 
a heuristic procedure. Since the filters are very lossy, it 
is desirable to increase e from the non-squeezed value, 
2ir x 25 Hz, so that the noise due to filter losses decreases 
and although the ideal minimum of g^ 03810138 increases, it 
is still buried by the noise due to filter losses. As we in- 
crease e from 2ir x 25 Hz, we search for the A and £ that 
minimize Sh at 600 Hz; we find that the sensitivity at 
600 Hz remains roughly the same, while the bandwidth 
increases. Trying to include as many sources as possible, 
we set e = 27r x 100 Hz for 5 dB squeezing and 27T x 60 Hz 
for 10 dB squeezing. The interferometer and filter param- 
eters used in these configurations are listed in Table. IVTT1 

As we see from Fig. the Harms et al. [two dark con- 
tinuous lines, one for 5dB squeezing the other for 10 dB 
squeezing] and the BC [two dark dashed lines] schemes 
are extremely close to each other. The peak sensitivities 
in the 5dB and 10 dB cases are chosen to be compara- 
ble to each other, while 10 dB squeezing gives a broader 
band. Although the FD techniques cannot increase the 
peak sensitivity much due to filter losses, they do in- 
crease the bandwidth of observation. This will allow the 
observation of multiple possible sources with a fixed con- 
figuration. For example, with the frequency and GW 
strengths estimates we used in Fig. [7] with 10 dB squeez- 
ing, we can detect simultaneously 7 sources near 600 Hz 
(including Sco X-l), while with 5dB squeezing we can 
detect 6 of them simultaneously (including Sco X-l). In 
Fig. HO we plot the increase in SNR by the squeezed 
schemes, as compared to the non-squeezed schemes, for 
LMXBs around the resonant frequency; both 5dB and 
10 dB squeezing are shown. In Table IVl columns 5-9 we 
list the sensitivies of these configurations. 

As in the case of NS binary inspirals, SR interfer- 
ometers with frequency-independent squeezing and read- 
out phase can also be optimized for the detection of 
LMXBs. However, squeezing combined with frequency- 
independent input-output optics cannot easily improve 
peak sensitivity and bandwidth at the same time for 
narrowband configurations. As a consequence, as we 
optimize the frequency-independent scheme with 5 dB 
squeezing, we obtain narrowband configurations that can 
detect at most 4 sources out of the group of 7 (includ- 
ing Sco X-l). [With 10 dB squeezing, when a similar 
optimization is done for frequency-independent schemes, 
one finds that a wideband interferometer with frequency- 
independent scheme 9 can detect all 7 sources — no nar- 
rowbanding is necessary, as we shall see in the next sec- 
tion.] 

Now we look at the interferometer performances if 



9 Corbitt, Mavalvala and Whitcomb !.'->— til are currently investigat- 
ing this optical configuration. 
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FIG. 9: Consequences of spin-frequency and emission-mechanism uncertainties on the detection of LMXBs (Z sources) using 
narrowband configurations. We plot the GW frequency and characteristic GW strength (for 20-day coherent integration), 
under approximations (al) (solid circles with center frequency around 600 Hz), (bl) (open circles, with center frequency around 
1200 Hz), (a2) (solid triangles, with center frequency around 400 Hz) and (b2) (open triangles, with center frequency around 
800 Hz), along with Harms et al. (equivalent to BC) noise curves (with 10 dB squeezing) tuned to those frequencies, with (solid 
curves) and without (dashed curves) internal losses of the interferometer. We assume an ITM power transmissivity of 0.033, 
SR-cavity round-trip loss of 1% and photodetection loss of 2%; 7 C = 840 kW and m = 40 kg. The rest of the parameters are 
fisted in Table IVTTT1 



assumptions other than (al) turn out to be true. In 
Fig. we show the predicted GW strengths from the 
Z sources under the four assumptions (obtained using 
Tables IVl and IVI|) : (al) (solid circles, with center fre- 
quency around 600 Hz), (bl) (open circles, with center 
frequency around 1200Hz), (a2) (solid triangles, with 
center frequency around 400 Hz) and (b2) (open trian- 
gles, with center frequency around 800 Hz). Given these 
hypothetical groups of sources, we tune squeezed-input 
SR interferometers (with Harms et al. or BC schemes, 
which are equivalent at high frequencies) to each of them: 
around 600 Hz, 1200 Hz, 400 Hz, and 8 00 Hz , with in- 
terferometer parameters listed in Table IVIIII and noise 
curves shown in Fig. [5] We remark that assumptions 
that yield lower /gw's tend to make the sources more 
detectable. In this study, we have also taken into ac- 



count interferometer losses, which has been neglected up 
till now. We assume the ITM power transmissivity to 
be T = 0.033, SR-cavity round-trip loss to be 1% (de- 
noted by Asr in Ref. |3l|) and photodetection loss to be 
2% (denoted by Apd in Ref. 31]). 10 These numbers are 
crude estimates; given the effects of interferometer losses 
suggested by Fig. [21 especially in higher frequencies [i.e., 
if assumptions (bl) or (b2) turns out to be true], more 
refined understanding of realistic interferometer losses, 



The 1% SR-cavity loss is the major interferometer loss, according 
to Ref. Li ill . In addition, we did not use the value T = 0.005 
in the baseline design of Advanced LIGO: assuming the same 
amount of loss per round trip inside the SR cavity, a much smaller 
T will make the effect of this loss much larger. 
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FIG. 10: Noise spectral densities of wideband configurations, 
with e = 2n x 600 Hz and A = 0, frequency-independent 
squeezing (a = n/2, £ = 0, light continuous curve), the 
Harms et al. scheme = 0, dark continuous curve) and the 
BC scheme (a = n/2, dark dashed curve) are used. Only the 
quantum-optical noise (taking filter losses into account) is in- 
cluded. Thermoelastic noises of spherical and MH mirrors 
are also shown, in dark and light dash-dot lines, respectively. 
The SQL is shown in dashed line. Possible GW signals from 
LMXBs [under the baseline assumption] and known radio pul- 
sars are also shown. 



as well as a more systematic study of interferometer pa- 
rameters will be crucial in fully understanding whether 
and how Advanced LIGO can detect these narrowband 
sources. 



C. Wideband configuration 

The so-called wideband configuration of SR interfer- 
ometers can be obtained setting A small and e rather 
high. These configurations can be used to detect a broad 
range of generic sources, including: coalescence of NS- 
NS binary, tidal disruption of NS by the BH companion, 
accreting NS's and radio pulsars. There are no specific 
criteria for the noise spectrum of the wideband configu- 
ration. For simplicity we set e = 2-7T x 600 Hz and A = 
(since this configuration is similar to the one by Corbitt 
and Mavalvala [23j], we denote it by CM). The various 
parameters used are summarized in Table llXl along with 
SNR achievable for NS binaries at 300 Mpc and sensitiv- 
ities at 600 Hz. Both 5dB and 10 dB squeezing are con- 
sidered, with 5 dB numbers quoted in parentheses, "[...]" . 
We plot the corresponding noise curves in Fig. 1101 

For frequencies higher than 200 Hz, the CM noise 
curves are always better than those with FD techniques. 
This is because, as observed by Corbitt and Mavalvala, 
at high frequencies, the optimal squeeze angle and detec- 
tion phase depend very mildly on the frequency. Thcrc- 



FIG. 11: Comparison of sensitivities to LMXB sources. We 
plot the noise curves of Harms et al. (or equivalently BC) 
scheme (continuous curves) and with frequency independent 
squeezing and readout scheme (the CM configuration, dashed 
curves), with 5dB (light curves) and 10 dB (dark curves) 
squeezing. We also plot the estimated frequencies and charac- 
teristic GW strengths of Sco X-l and several other Z sources 
[under the baseline assumption (al)]. 



fore, the FD schemes, having additional filter losses, give 
worse performances. At high frequencies, the wideband 
schemes give a sensitivity of 3.2 (10 dB squeezing) or 1.8 
(5 dB squeezing) times better (in amplitude) than the 
wideband configuration without squeezing. With 10 dB 
squeezing, the wideband configurations can detect known 
pulsars at lOkpc with e > 10~ 7 if /gw ~ 420 Hz, with 
e > 3 x 10~ 8 if f GW > 1 kHz. [With 5 dB squeezing, the 
minimum detectable e will be 1.8 times larger than the 
10 dB value.] However, if we also require good sensitiv- 
ities below 200 Hz, then the FD wideband schemes are 
preferable to the CM configuration . 

In addition, in the 10 dB squeezing case, when spheri- 
cal mirrors are used, the SNR for binaries are all above 
96% the optimal values obtained in the broadband case 
(see Table Hvj) . However, for Mexican-Hat mirrors, the 
SNR is less optimal, equal to 83% (no filters), 91% 
(Harms et al.) and 93% (BC) the optimal values (of 
the same scheme). [See Table llVl ] These can be under- 
stood by going back to Sec. IV Bl and observing in (the left 
panel of) Fig. [5] that for spherical mirrors, the optimal 
noise curves are very wideband. 

It is also interesting to note that, with 10 dB squeez- 
ing, the sensitivities of wideband configurations around 
600 Hz, are only slightly worse, ~ 10% in amplitude, 
than the narrowband configurations. As a consequence, 
with 10 dB squeezing, the wideband configurations, even 
without FD techniques, can detect the same groups of 
LMXBs discussed in the last section (see Fig. HOfl . How- 
ever, it should be noted that, if 10 dB squeezing is not 
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Interferometer Configuration 


Filter I 


Filter II 


Performance 10 dB [5 dB] 


e X 

Scheme a C 

2ttHz 2ttHz 


Q ™ T] 
27r Hz (ppm) 


°" a if 
27rHz (ppm) 


V^/(10~ 24 /\/H5) NS-NS at 300 Mpc 
at 600 Hz Sphcrical/MH 


No Filters 600 tt/2 
Harms et al. 600 FD 
BC 600 tt/2 FD 


2.8 - 600.0i 1508 
2.8 - 600. Oi 1508 


-41.3- 41. Oi 103 
-41.3- 41. Oi 103 


0.99 [1.77] 6.47/8.70 [5.69/7.60] 

1.12 [1.82] 7.00/10.65 [6.01/8.76] 

1.13 [1.85] 6.68/9.68 [5.79/8.09] 



TABLE IX: Parameters of wideband configurations, including: squeezing with frequency independent input-output optics, 
the Harms et al. scheme and the BC scheme. Since there is no detuning in these SR interferometers, they are equivalent to 
conventional interferometers, which KLMTV studied. We also give the noise spectral densities at 600 Hz, and the SNR for 
NS-NS binaries at 300 Mpc distance. Both 5dB and 10 dB squeezing are considered, with 5dB numbers quoted in parentheses, 
"[...]". In cases with 10 dB squeezing, the sensitivties at 600 Hz is only slightly worse than the narrowband configurations (see 
Table EH- 
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FIG. 12: Consequences of spin- frequency and emission- 
mechanism uncertainties on the detection of LMXBs (Z 
sources) using the wideband configuration. We plot the GW 
frequency and characteristic GW strength (for 20-day coher- 
ent integration), under approximations (al) (solid circles), 
(bl) (open circles), (a2) (solid triangles) and (b2) (open trian- 
gles), along with CM noise curves with 5dB (dashed curves) 
and 10 dB (continuous curves) squeezing, with (dark curves) 
and without (light curves) interferometer losses. 



achievable, then one cannot detect these sources with 
the wideband configuration. For example, 5 dB squeez- 
ing will barely allow one or two more LMXBs than Sco 
X-l to be detected. The narrowband configuration (with 
FD input-output schemes), by contrast, will only miss 
one source in the group of 7. In Fig. ^2 we compare the 
sensitivities of narrowband FD schemes and wideband 
frequency independent schemes to LMXB sources, with 
5dB and 10 dB squeezings. 

Finally, by taking into account all other assumptions 
on spin frequency and GW emission mechanism, we plot 
in Fig. U3 the predictions of (al), (bl), (a2) and (b2), 
along with CM noise curves with 5 dB (dashed curve) and 
10 dB squeezing (continuous curve), with (dark curves) 
and without (light curves) interferometer losses included. 



VI. THIRD-GENERATION 
INTERFEROMETERS 

We now assume that on time scales of third-generation 
GW interferometers (around 2012), thermal noise of mir- 
rors will be reduced by a large factor, for example by us- 
ing cryogenic techniques, and we can take full advantage 
of the improvements in quantum noise obtained by FD 
input-output techniques. In addition, we assume that 
long filters can be fit into the existing vacuum tubes 
(which house the arm cavities) of the LIGO facility and 
made 4 km long, so that optical losses will be significantly 
lowered (see Table ILTiTl . 



As discussed in Sec. IIVDI the 
BC scheme is nearly optimal for frequencies lower than 
the optical resonance [see Fig. Eland Eq. thus in the 

following we shall restrict our analysis to the BC scheme. 
However, before showing the performances, we want to 
discuss the limitations of the so-called short-cavity ap- 
proximation, so far used in the literature to describe 
kilometer-scale filter cavities 20, 261. 



A. Breakdown of short-cavity approximation 

Up till now in this paper, we have been using the short- 
cavity approximation, which imposes that flL/c <C 1. 
[Note that when refered to the interferometer, L is the 
arm length, ft is the GW sideband frequency or the op- 
tical resonant frequency —A — ie; when refered to filter 
cavities, L is the filter length, is the GW sideband fre- 
quency or the filter resonant frequency O reB .] As we saw 
in Sees. Inland IIIII the short-cavity approximations, ap- 
plied to SR interferometers and KLMTV filters, simplify 
significantly their input-output relations [see Eqs. J3J) - 
(JSJ), (|38fl and (|42|) ]. allowing a straightforward determi- 
nation of filter parameters in the Harms et al. and BC 
schemes via characteristic equations [Eqs. (|7TJ) and l|74|)]. 

On the contrary, without this approximation (i.e., 
when cavity lengths are too long for this approxima- 
tion to work), the filter parameters cannot be deter- 
mined easily — it is not even clear whether the opti- 
mal/suboptimal frequency dependence required by (the 
exact input-output relation of) SR interferometers can at 
all be realized by (those of) KLMTV filters. 

Since we have derived the exact input-output relation 
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FIG. 13: To show the break-down of the short-cavity ap- 
proximation we plot several noise curves for the KLMTV 
(conventional) squeezed-variational interferometer [2(J with 
I c = 840 kW, m = 40 kg and e" 2r = 0.1, fixing a bandwidth 
of e = 2tv x 75 Hz. The dark continuous curve refers to the 
nominal filter-cavity length Lf = 4 km, and the round-trip 
filter loss T e = 20ppm; while the short-filter approximation 
(Lf = 0.1 m) predicts the dark dashed curve. For comparison, 
we also plot noise curves of configuration with Lf = 400 m, 
T e = 2ppm (dark dash-dot curve), along with those of the 
lossless optical configurations with Lf = 4 km (light continu- 
ous curve), Lf = 400 m (light dash-dot curve) and with short- 
filter approximation (Lf = 0.1 m) (light dashed curve). 



of the filters [Eqs. as well as (partiall y 11 ) that 

of the interferometer [Eqs. (99)-(104) of Ref. |2£|], we 
can investigate the range of validity of the short-cavity 
approximations. 

Let us start with conventional interferometers. As we 
have checked in this case, the short-arm approximation is 
still quite accurate, in the sense that, for a given readout 
scheme (i.e., a given set of input or output filters), using 
exact and short-arm-approximated input-output relation 
do not give very different results. Yet, the short-filter 
approximation seems to lose accuracy at low frequencies. 
We study this effect in Fig. ^| by plotting several noise 
curves for squeezed-variational conventional interferom- 
eters HJ with I c = 840 kW, m = 40 kg, e - 2r = 0.1, 
and e = 2?t x 75 Hz, using the exact interferometer input- 
output relation. In doing so, we use filters with band- 
widths and resonant frequencies obtained from the short- 
filter approximation, but with different actual lengths 
and losses. In the figure, we show the noise curve for 
filters with Lf = 4 km and T e = 20 ppm in dark contin- 



11 In Ref. l3dl we treated exactly the propagation of light inside 
the interferometers, approximated the radition-pressure-induced 
motion of the ITM as being equal to that ot the ETM. 



FIG. 14: To investigate the short-filter and short-arm approx- 
imation in SR squeezed-variational interferometers we plot 
several noise curves fixing I c = 840 kW, m = 40 kg, e~ 2r = 
0.1, e = 2n x 80 Hz and A = 2n x 200 Hz. In particular we 
show the noise curve obtained with exact input-output rela- 
tion and 4 km filter cavities (dark continuous curve), and 4 m 
filter (dark dashed curve) ; the noise curves obtained with first- 
order expanded input-output relation and 4 km filter (light 
continuous curve) and 4 m filter (light dashed curve) . We use 
T e = 20 ppm for the 4 km configurations, and T e = 0.02 ppm 
for 4 m configurations, such that the overall loss factor re- 
mains the same. The disagreements between curves with the 
same pattern (continuous or dashed) but different color (dark 
or light) is due to the inaccuracy of short-filter approximation; 
while the disagreements between curves with the same color 
but different pattern is due to the inaccuracy of short-arm 
approximation. 



uous curve, and also lossy filters with decreasing length 
but the same T e /Lf ratio: Lf = 400 m in dark dotted 
curve and Lf = 0.1m (to simulate short-filter limit) in 
dark dashed curve. The noise spectrum improves as the 
filter length decreases. [In fact, since in this case the 
short-arm approximation is accurate, short filters must 
give the optimal performance.] In Fig. 1131 we also show 
noise curves for lossless configurations with Lf = 4 km in 
light continuous curve, Lf — 400 m in light dotted curve 
and Lf = 0.1m (to simulate short-filter limit) in light 
dashed curve. The reason for such dramatic noise in- 
crease at low frequencies can be attributed to the strong 
ponderomotive squeezing generated by conventional in- 
terferometers at these frequencies (note that q — > +oo 
as CI — > 0, see left panel of Fig. The stronger the 
squeezing, the higher the accuracy requirement on the 
FD readout phase; yet the accuracy of short-filter ap- 
proximation does not increase indefinitely when CI — ► 0. 

By contrast, as we have checked, the short-cavity ap- 
proximations still apply very well to squeezed-input con- 
ventional interferometer which at low frequencies does 
not have as good an ideal sensitivity as the squeezed- 
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variational conventional interferometer. 

In Fig. ^] we investigate the short-arm and short-filter 
approximations for SR squeezed-variational interferom- 
eters (the BC scheme) with I c = 840 kW, m = 40 kg, 
e~ 2r = 0.1, e = 2tt x 80 Hz and A = 2tt x 200 Hz. In 
this case, both the short-filter and short-arm approxima- 
tions introduce some inaccuracies, but they are by far not 
as significant as in the squeezed-variational conventional 
interferometers. In particular, in Fig. 1141 we plot noise 
curves obtained using exact interferometer input-output 
relation, with 4 km (dark continuous curve) and 4 m fil- 
ters (dark dashed curve) , and noise curves obtained using 
short-arm-approximated interferometer input-output re- 
lation, with 4 km (light continuous curve) and 4 m (light 
dashed curve) filters. We fix T e = 20ppm for 4km con- 
figurations, and 0.02 ppm for 4 m configurations, keeping 
the same overall loss factor. [Filter resonant frequencies 
and bandwidths arc still obtained from the characteristic 
equation (|74() . which in turn has been derived based on 
both short-arm and short-filter approximations]. Noise 
curves with the same color (light or dark) use the same 
interferometer input-output relation, so the difference be- 
tween them reflects the inaccuracy of the short-filter ap- 
proximation; those with the same pattern (continuous 
or dash) share the same filter input-output relation, so 
their difference reflects the inaccuracy of the short-arm 
approximation. We conclude that the errors arising from 
the short-arm and short-filter approximations somewhat 
cancel each other, making the exact noise curve differ 
only slighly from the curve with both short-arm and 
short-filter approximations applied. The mild noise in- 
crease around the optical-spring resonance in this case 
can also be understood from the ponderomotive squeez- 
ing factor. As we see from the left panel of Fig. (the 
dashed curve represents a similar configuration), pon- 
deromotive squeezing is the strongest near this resonance, 
yet even here the squeeze factor is still small compared 
to that of the conventional interferometer at lower fre- 
quencies. 

We now discuss the short-cavity approximation in 
squeezed-variational speed meters [2^, |27j. We con- 
sider the configuration with £! s = 2ir x 95.3 Hz (the 
"sloshing frequency", as denoted by Q in Ref. |2l|) and 
S = 2n x 100 Hz (bandwidth), assuming e~ 2r = 0.1. We 
include optical losses as done in Ref. |21|. As for con- 
ventional squeezed-variational interferometers, the short- 
arm approximation is rather accurate here. [This is true 
if the enhanced formula (i.e., expanded to next-to- leading 
order in CIL/c) for the quantity k is used, see footnote 5 of 
Ref. |21| .] However, the short-filter approximation is not 
accurate enough, if we increase the optical power further 
from Advanced LIGO value. In Fig. ED we plot four noise 
curves with I c = 840 kW (dark curves) and 2 MW (light 
curves), and filter lengths 4000m (and T e = 20 ppm, 
continuous curves) and 4 m (and T e — 0.02 ppm, dashed 
curves). [Again, resonant frequencies and bandwidths of 
the filters are obtained in the same way as in Ref. |2lj . 
based on short-arm and short-filter approximations.] As 
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FIG. 15: We plot the noise curves of squeezed-variational 
speed meters with Q s = 2ir x 95.3 Hz, 8 — 2n x 100 Hz, and 
e~ 2r = 0.1, assuming I c = 840 kW (dark curves) and I c = 
2MW (light curves), and Lj — 4000 m (continuous curves) 
and Lf = 4m (dashed curves). The round-trip losses are 
T e = 20 ppm for 4000 m filters and T e = 0.02 ppm for 4 m 
filters. The optical losses are included following Ref. |2l||. 



we sec, a filter length of 4000 m increases the noise signif- 
icantly as I c becomes on the order of 2 MW. The increase 
is rather constant (and now as dramatic as in KLMTV 
squeezed-variational conventional interferometers) at low 
frequencies, because speed meters have a constant pon- 
deromotive squeezing factor at low frequencies |2lj . 

We notice that in all the above cases where the short- 
cavity approximations break down, using filter param- 
eters obtained from the characteristic equations (as we 
have done above), which are derived assuming those ap- 
proximations, can no longer be optimal. Instead, one 
must optimize filter parameters numerically using exact 
filter and interferometer input-output relations. We do 
not have quantitative results yet on how much sensitiv- 
ity can be gained by this re-optimization, but it does 
not seem likely that the sensitivity can reach the opti- 
mal level (i.e., having the FD rotation from the filters 
matching exactly the interferometer's requirement). 



B. Performances of SR squeezed-variational 
interferometers 



Using exact filter and interferometer input-output re- 
lations (i.e., without applying short-cavity approxima- 
tions), and assuming that 4 km filters will be used in 
third-generation interferometers, we compare in Fig. 1161 
the noise spectral densities of conventional squeezed- 
variational interferometers, SR squeezed-variational in- 
terferometers (BC scheme), and the squeezed-variational 
and -input speed meters. The BC scheme (which re- 
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FIG. 16: Noise spectral densities ol SR interferometers with 
BC scheme (SR+BC, dark continuous line), conventional 
squeezed-variational interferometer (KLMTV, dark dashed 
line), the Purdue-Chen speed- meter interferometer with or- 
dinary homodyne detection (light continuous line) and FD 
homodyne detection (light dashed curve). We assume I c = 
840 kW, m = 40 kg and e~ 2r = 0.1; the SR interferome- 
ter with BC scheme has A = 2ty x 200 Hz, e = 2n x 80 Hz 
and a — arctan(e/A); the speed meter with ordinary homo- 
dyne detection has Q s = 173.2 Hz and S — 2tt x 200 Hz, 
while the speed meter with FD homodyne detection has 
Sl s = 27T x 95.3 Hz and S = 2ty x 100 Hz. Optical filters are 
assumed to be 4 km long, with 20ppm round-trip loss. 



quires two additional km-scale cavities) has better sensi- 
tivity than the conventional squeezed-variational scheme 
(which also requires two additional km-scale cavities) for 
all frequencies below ^350 Hz. It has also better per- 
formances than the squeezed-input speed meter (which 
requires one additional km-scale cavity) for all frequen- 
cies above ~40Hz. The BC scheme has comparable (or 
slightly better) sensitivities with respect to the squeezed- 
variational speed meter (which requires three additional 
km-scale cavities) 12 for frequencies between ~ 50 Hz and 
- 300 Hz. 



VII. CONCLUSIONS 

In this paper, we generalized the study of KLMTV [2(j 
on FD input-output optics to SR interferometers, and 
discussed possible applications to second- and third- 
generation GW interferometers. 



We do not discuss the Sagnac interferometer which is also a speed 
meter without adding any km-scale cavities 1271 . A Sagnac in- 
terferometer can achieve sensitivities equivalent to the Michelson 
Purdue-Chen speed meters, and its squeezed-variational version 
requires only two additional km cavities. 



In the first part of the paper (Sees. ITT1- IIV|) . we stud- 
ied the quantum optical properties of SR interferome- 
ters and FD input-output schemes. We wrote the input- 
output relations of SR interferometers as a product of 
ponderomotive squeezing and quadrature rotations, de- 
riving explicit formulas for the intrinsic rotation angle 
and squeeze factor [see Eqs. <|20|) and (|21|) ]. and inves- 
tigating their features for several optical configurations. 
We found that ponderomotive squeezing becomes very 
weak in SR interferometers for frequencies higher than 
~ 300 Hz, regardless of the optical configuration [see 
Eq. i|25|) ]. Then, we built and analyzed the performances 
of the input-output scheme which combine FD homo- 
dyne detection (via KLMTV filters) with ordinary input 
squeezed vacuum (BC scheme), and compared it to the 
recent FD scheme proposed by Harms et al. [25|. In the 
low-power limit (which also describes the high-frequency 
band of Advanced LIGO) we worked out the fully opti- 
mal input-output scheme [see Eq. I|83l) ]. In the general 
case, we derived simple analytical formulas for the fully 
optimal noise spectrum [Eq. I|90|l ] and the optimal input 
squeeze angle [Eq. (|89|) ]. and found that at low frequen- 
cies, the BC scheme can approximate the fully optimal 
noise curve very well [see Eq. (I93|) ]. providing better per- 
fomances than the Harms et al. scheme. These results for 
SR interferometers are quite similar to the conventional 
interferometer case, in which as shown by KLMTV, a 
frequency independent squeezed vacuum is already fully 
optimal (with FD readout) , yet a frequency independent 
readout cannot give as good a sensitivity (even with FD 
squeezing). [The BC and Harms et al. schemes gener- 
alize to SR interferometers the squeezed-variational and 
squeezed-input schemes introduced by KLMTV for con- 
ventional interferometers.] 

In the second part of the paper (Sec.EJ, assuming that 
squeezed vacuum in the GW band would become avail- 
able during the operation of Advanced LIGO, we evalu- 
ated the improvement in astrophysical sensitivity to spe- 
cific sources achievable by these FD schemes, under the 
facility limitation that the filters cannot be longer than 
30 meters. It is important to note that, as has been real- 
ized by Corbitt and Mavalvala [23[, for nearly tuned SR 
interferometers with a large bandwidth (wideband con- 
figuration), the optimal input-output scheme is nearly 
frequency independent at high frequencies. So, in this 
case it is possible to use squeezing optimally without in- 
troducing FD techniques. The Corbitt-Mavalvala (CM) 
wideband configuration can be used to detect simulta- 
neously various types of sources in the high-frequency 
band, e.g., NS-NS merger, tidal disruption in NS-BH sys- 
tems, or GWs from known radio pulsars. In addition, if 
10 dB squeezing can be realized, this wideband configu- 
ration can already detect a group of 7 LMXBs (including 
Sco X-l) around 600 Hz. 

However, for specific sources with known spectral fea- 
tures, it is more convenient to use optical configurations 
which are not wideband. In this case, the FD techniques 
can provide more flexibility and somewhat better sensi- 
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tivity, despite significant optical losses due to short filters. 
The Harms et al. scheme is shown to provide a better sen- 
sitivity than the BC scheme in general in the Advanced 
LIGO era, due to the BC scheme's higher susceptibility 
to losses at low frequencies and the filter-length limita- 
tion to 30 meters. For NS-NS inspirals, assuming 10 dB 
squeezing, without using FD filters, one can improve the 
event rate by 89% (spherical mirror) or 42% (MH mirror) 
with respect to the non-squeezed case; further improve- 
ments of 20% (spherical mirror) and 42% (MH mirror) 
can be obtained by the Harms et al. scheme. For LMXBs, 
using 5 dB squeezing, without using FD techniques, the 
broadband configuration can only detect 3 sources si- 
multaneously. By adjusting the SR parameters, we find 
that a frequency independent input-output scheme can- 
not detect more than 4 sources, since in this case we 
cannot gain sensitivity and bandwidth at the same time. 
The Harms et al. and BC schemes, instead can allow 
the detection of 6 sources simultaneously, by opening up 
the bandwidth with FD filters (although the peak sensi- 
tivity cannot be improved much due to significant filter 
losses). With 10 dB squeezing, the Harms et al. and BC 
schemes can open up the bandwidth further, including a 
group of 7 sources near 600 Hz. [However, with this level 
of squeezing, the broadband configuration can also de- 
tect the same sources, though with slightly less sensitiv- 
ity.] It is important to mention that, the above detailed 
results about LMXBs are obtained by assuming (base- 
line assumption) that the Z sources, which are the most 
promising GW sources among all the LMXBs, are spin- 
ning at the QPO difference frequency (f s — /<j ~ 300 Hz) 
and that the GW is dominantly mass-quadrupole emis- 
sion (/gw = 2/ s 600 Hz). Recent results [45|, however, 
may suggest that the spin frequencies could be twice the 
QPO difference frequency and therefore lie around 600 Hz 
instead. In addition, it is not clear yet whether the domi- 
nant GW emission is mass quadrupole (e.g., due to defor- 
mation in the crust) or current quadrupole (e.g., due to 
r-mode). In Figs. 151 and !12l we have briefly explored the 
sensitivities of narrowband and wideband configurations 
to these alternative scenarios. Moreover, we also realized 
that optical losses inside the SR interferometer become 
rather crucial in deciding whether these sources will be 
detectable, especially if the predicted waves are at higher 
frequencies. A more careful study of optical losses will be 
reported in a forthcoming paper 50]. It should also be 
remembered that we have been relying on the character- 
istic GW strength j3|| of LMXBs to provide a very rough 
criterion of detectability — in order to make a straight- 
forward yet specific comparison between different noise 
curves. As a consequence, the true detectability of these 
sources by configurations studied in this paper should be 
refined by a more rigorous study. 

In the third part of the paper (Sec. IVII) . we investi- 
gated the performances of squeezed SR interferometer 
with FD output using two 4-km KLMTV filter cavities. 
We found [see Fig. that SR interferometers with in- 
put squeezing and sub-optimal FD readout scheme (the 



BC scheme) can have competing sensitivities to existing 
proposals for third-generation interferometers [2^, [2^, es- 
pecially in the middle frequency band of 50 - 350 Hz (see 
Fig. 1161) . We also discuss the limitation of the short-arm 
and short-filter approximations adopted by most of the 
past works [53, |5g . Should at least one of these approx- 
imations break down, the optimal (or suboptimal) filter 
parameters provided by the characteristic equation l|39|l 
would not give the required FD rotation — which, in 
fact, may not even be realizable by any sequence of de- 
tuned FP cavities. Indeed, we found that for squeezed- 
variational conventional interferometers and speed me- 
ters (with high power), the approximation breaks down 
at low frequencies if applied to km-scale filters (see Fig. 1131 
and Fig. 1150 . which we attribute to the high level of 
ponderomotive squeezing and therefore more stringent 
requirement rotation accuracy in these configurations; 
while for SR squeezed-variational interferometers (BC 
scheme), the error introduced by the approximations is 
rather mild (see Fig. I14f> . 
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APPENDIX A: PROOF THAT THE FULLY 
OPTIMAL SCHEME SATISFIES THE 
SUB-OPTIMAL CONDITION 

Suppose £ is the readout phase. As seen in Sec. 
IV, with the sub-optimal a given by Eq. I|67|l . only the 
squeezed quadrature enters the detected quadrature. We 
now show that given £, this a is also the optimal squeeze 
quadrature in the sense that it minimizes the output 
noise. Since when we vary a alone, the signal strength in 
the output quadrature remains constant, we only need 
to minimize the noise in b^. We write Eq. (|77|l schemat- 
ically as 

Ab ( = {A 1 A 2 )R(a)( 6 ^|) , (Af) 

with (Ai A2) the product of the readout and the inter- 
ferometer parts, Ai t 2 € 5ft. The noise spectrum depends 
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on a as 



S h oc ( At A 2 ) R(a) 



-2r 



Minimizing S/j then requires 



(A2) 



(A 1 An) R(a) J J - , (A3) 
which is equivalent to Eq. (|78|l and hence to Eq. (|67[) . 



APPENDIX B: ON THE DETECTABILITY OF 
NARROWBAND SOURCES 

In this appendix, we briefly review some subtleties 
that are not taken into account in our discussion of nar- 
rowband sources. We restrict the analysis to the mass- 
quadrupole radiation mechanism. 



the SNR over different source and detector orientations. 
If we want to understand the detectability of a particu- 
lar source and extract an upper limit, we should consider 
unfavorable geometries. In this paper the predictions for 
LMXBs and isolated pulsars have been obtained averag- 
ing the SNR. As said, this is not appropriate for evalu- 
ating the detectability of individual sources. In what fol- 
lows we shall briefly review the average procedure, and 
comment on what might be done in order to extract up- 
per limits from a specific source. 

Let us first consider the variation or uncertainty in 
f+,x- The most straightforward ansatz for taking this 
into account is to use the r.m.s. average of SNR over the 
entire sky — as viewed by the detector. The ansatz gives: 



(SNR 2 



Met 



V2(F_ 




Characteristic amplitude of monochromatic 
gravitational waves: h c 



To obtain the RHS in the above equation we use (F]_) = 
(F 2 ) = 1/5, and define [as done in Eq. (51) of Ref. |5lj]: 



Let us consider a monochromatic source emitting GWs 
at (angular) frequency fig : 



h(t) = F + h + cosftot- 
h + = h (l + cos 2 i), 
h x = 2ho cosi , 



F x h x sin Slot , 



(Bl) 
(B2) 
(B3) 



where we denote with F + , F x the antenna patterns [5l| 
and with i the angle that the line-of-sight forms with the 
spin direction of the neutron star. The quantity ho is an 
intrinsic GW amplitude depending on the ellipticity of 
the isolated pulsar, or on the X-ray flux emitted by the 
LMXB (through the balance between accretion torque 
and GW radiation-reaction torque), as well as on the 
distance of the source. Suppose the signal is observed in 
the time interval —To/2 < t < To/2, then the signal-to- 
noise ratio (SNR) using optimal matched filtering is: 



SNR 



+ 00 



drt2\h(n)\ 2 



2tt Sh(fi) 



2Tn 



i r p2 



S h (Qo) 



h* + +F*h% 



(B4) 



Note that Eq. (|B4j| differs from Eq. (29) of Ref. |E3 - 
the latter is wrong by a factor of 2 [52j . As we shall see in 
the next section, when statistical issues are considered, 
signals above a certain threshold SNR will be detectable. 

In reality i*+, x vary due to earth's motion, or are 
unknown for some sources; the inclination angle i, al- 
though may stay constant, could also be unknown, even 
for known sources, e.g., LMXBs. As a consequence, with 
a fixed ho, the SNR achievable can be different. If we 
are interested in expected event rates, we should average 



h n (ft ) 



is h (n )/T 



2{F%) 



(B6) 



Now let us consider the dependence of y (SNR 2 ) dot on 
i. There are two plausible averaging prescriptions. The 
first, which is the easiest, and most appropriate for a 
known source at a fixed distance, averages SNR 2 uni- 
formly over source angles, as 



(SNR 2 )dot&src = T 




It is then natural to define as characteristic amplitude: 



ho ~ 1.26 ho . 



(B8) 



Another way of averaging was proposed by Thorne in 
Ref. [SJ, which has the property that event rate oc 
(SNR) 3 , if we assume uniform distribution of sources in 
the universe, 



(SNR) 



dot & sre 



c?f2 s 



4tt 



[(SNR 2 ) dct ] 



1 

h n 



47T 



h% + ^ 



(B9) 
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The integral l|B9(l cannot be performed analytically, so 
Thorne introduced a kludge factor 



dfl SIC fh\ + h 2 ^ 2 



4tt 
yielding 



dft aTC h\ + h\ 



4tt 



[(SNR)dct&src] k i udgc = ^ y ^ 




(BIO) 



(Bll) 



The above expression originated the following definition 
for the characteristic strength: 



^cj kludge 



1.46 h - 



(B12) 



The kludged characteristic strength [/i c ]kiudge nas been 
used by many authors, including us in this paper. In 
particular, [/i c ] kludge of a pulsar at distance r with ellip- 
ticity e and frequency / can be obtained from Eq. (3.6) 
of Ref. [3^, and [/i c ] kludge from LMXBs with mass- 
quadrupole emission mechanisms that balances the ac- 
cretion torque can be obtained from Eq. (4) of Ref. • 
However, the kludge factor \/4/3 is not accurate. A 
simple numerical calculation gives: 



3 1 

2 1 h2 \ 2 



dfl src ( hi. + h 

47T 



dn src h\ + h% 



4tt 



sini [(1 + cos 2 if + (2 cosi) 2 ] 3/2 di 



1 

2 Jo 
1.047. 



sini [(1 + cos 2 i) 2 + (2cosi) 2 ] di 



h c = 1.32 ho . 



(B13) 



(B14) 



On the other hand, if we are interested in setting upper 
limits, we should use an i that has the lowest possible 
(SNR)dct; this implies i — ir/2 and 



min {(SNR)det} = , h Uh = J^h = 0.7Q7/i . 

(B15) 

For a known source with constant ho but uncertain 
orientation (uniformly distributed cosi), we can also 
ask for the probability that (SNR)dct exceed h a /h n , 
[h c ] kludge /h n and h c /h n . The answers are 41%, 29% and 
37%, respectively. 

To summarize, we have managed to write the SNR 
for a given source or a given set of sources (with fixed 



intrinsic amplitude ho, unknown i) in the form of 

h c . 



SNR= -± 
h n 



(B16) 



where h c is the characteristic amplitude — with four 
different relations to h , {B8j|, {Bl]|, (|B"T4|) and jHTSJl, 
yielding SNRs that are either averaged in different ways 
over different i's, or taken as the minimum. [The quan- 
tify h n is defined in Eq. (|B6(I in terms of Sh and inte- 
gration time.] Given the characteristic amplitude based 
on a particular prescription, values based on other pre- 
scriptions can be obtained from Eqs. (|B8|I . (|B12() . (|B14() 
and l|B15(l using the fact that ho is the same in all of 
them. For example, given [/i c ]kiudge (which is used in 
this paper), we have 



h a = 0.866 [/ijkiudgc : 

h c = 0.907 [/ijkludge : 
/lUL = 0.484 [/ijkludge • 



(B17) 
(B18) 
(B19) 



Note that the more conservative /iul is a factor of 
smaller than [h c ] kludge- 



2. Factors that determine the detection threshold: 
from h c to Sh c 

With data analysis methods and desired statistical 
confidence, a threshold (minimum) SNR can be obtained; 
hence from Eq. (|B16|) . for a certain h c , a maximum h n , 
and thus a maximum Sh, or Sh a can be obtained. [Here 
h c should be specified from an intrinsic GW amplitude 
ho, through a h c -ho relation, like one of Eqs. i|B8(l . I|B12(I . 
(|BTIj) and JH33.] 

Refs. j33, y3 introduces the canonical sensitivity 
/13/yr = 4.2 \/S n (f) x 10 _7 Hz, which is the characteris- 
tic amplitude of the weakest source detectable with 99% 
confidence level (i.e., 1% false alarm) in a coherent search 
of 10 7 s of data, if the frequency and phase evolution of 
the source is known [see Eq. (1.4) of Ref. |38j]. This read- 
ily gives the Sh c for known pulsars. However, it should be 
noted that Sh c obtained using this h 3 / yr only guarantees 
that the expectation value (or average) of the detection 
statistic be higher than the detection threshold |3g. and 
gives a high false-dismissal rate of about 50% (a, |53j . In 
other words, even if the noise curve touches Sh c for a 
particular source, there is still around 50% chance this 
source will not make the detection threshold. 

For LMXBs, frequency and phase evolution of the GW 
due to orbital motion is unknown and one must build 
an appropriate bank of templates to search for these 
parameters, resulting in a threshold higher than h 3 / yl ; 
in addition, variations in accretion rate, which induces 
"random walks" in the spin, and hence in the GW fre- 
quency, further complicates the data analysis procedure, 
increasing the threshold further. Brady and Creighton 
studied these issues, and devised a two-step hierarchi- 
cal scheme for detecting such signals [3^. They use 
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the relative sensitivity re i to measure the increase in 
the threshold for the characteristic amplitude: = 
/i3/yr/©rei- As a consequence, in our notation, we 
have y/S hc (f) = h c 6 rc i/(4.2 VlO- 7 Hz). Brady and 
Creighton have shown that, for Sco X-l, with realistic 



computational power, 8 rc i = 0.41 [Sec. VIIC of Ref. |39j]. 
This yields a value of S/ lc comparable to that of a coher- 
ent integration of 20 days. In the paper we use this pre- 
scription for all LMXBs (except for SAX J1808. 4-3658) , 
and denote this characteristic strength by day . 
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